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1 - Introduction 

The aim of this paper is to derive an error estimate for approximate solutions of the compressible 
barotropic Navier-Stokes equations obtained by a discretisation scheme. These equations are posed on 
the time-space domain Qt = (0, T) x 0, where 0 is a bounded polyhedral domain of M. d , d = 2,3 and 
T > 0, and read: 


dtQ + div(fm) = 0, (1.1a) 

dt(gu) + div(fm <g) it) — pAu — (p + A)V div-u. + W x p(g) = 0, (1.1b) 

supplemented with the initial conditions 


£>(0,x) = Qo(x), gu(0,x) = q 0 u 0 , 


( 1 . 2 ) 


where £>o and uq are given functions from 0 to R+ and M d respectively, and boundary conditions 


u \{o,T)xdn — 0 - ( 1 - 3 ) 

In the above equations, the unknown functions are the scalar density field g(t, x) > 0 and vector velocity 
field u = (tti,.. ., Ud)(t, x), where t € (0, T) denotes the time and x G 11 is the space variable. The 
viscosity coefficients p and A are such that 

p > 0, A + —p T 0. (1.4) 

The pressure p is a given by an equation of state, that is a function of density which satisfies 

p G C([0, oo)) n C' 1 (0, oo), p( 0) = 0, p'{g ) > 0. (1.5) 


In addition to (USD, in the error analysis, we shall need to prescribe the asymptotic behavior of the 
pressure at large densities 


lim 


p'{q) 


irf = Poo > 0 with some 7 > 1; 


e-> 00 gP~ 

furthermore, if 7 < 2 in (USD , we need the additional condition (for small densities): 

p' (g) 

lim inf —— = po > 0 with some a < 0. 


0 g' 


i«+i 


( 1 . 6 ) 


(1.7) 


The main underlying idea of this paper is to derive the error estimates for approximate solutions 
of problem (ED-Q obtained by time and space discretization by using the discrete version of the 
relative energy method introduced on the continuous level in nanang. In spite of the fact that the 
relative energy method looks at the first glance pretty much similar to the widely used relative entropy 
method (and both approaches translate the same thermodynamic stability conditions), they are very 
different in appearance and formulation and may provide different results. The notions of relative en¬ 
tropy and relative entropy inequality were first introduced by Dafermos [7J in the context of systems 
of conservation laws and in particular for the compressible Euler equations. The relative energy func¬ 
tional was suggested and successfully used for the investigation of the stability of weak solutions to 
the equations of viscous compressible and heat conducting fluids in [13]. In contrast with the relative 
entropy of Dafermos, for the viscous and heat conducting fluids, the relative energy approach is able 
to provide the structural stability of weak solutions, while the relative entropy approach fails in this case. 


Both functionals coincide 


(modulo a change of variables) in the case of (viscous) compressible flows 
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in the barotropic regime. The relative energy functional and the intrinsic version of the relative energy 
inequality have been recently employed to obtain several stability results for the weak solutions to these 
equations, including the weak strong uniqueness principle, see H2JEE5]. Note that particular versions of 
the relative entropy inequality with particular specific test functions had been previously derived in the 
context of low Mach number limits, see e.g. [[521134] . 

The discrete version of the Dafermos relative entropy was employed in the non viscous case to de¬ 
rive an error estimate for the numerical approximation to a hyperbolic system of conservation laws and, 
in particular, to the compressible Euler equations [4]. In this latter paper, the authors assume an L°° 
bound for the discrete solution, which is uniform with respect to the size of the space and time disretiza- 
tion (usually called stability hypothesis), that is not provided by the discrete equations. The same 
method with the same severe hypotheses have been used in [35j to treat the compressible Navier-Stokes 
equations. The error analysis in the present paper relies on the theoretical background introduced in |T2] 
and yields an unconditional result, in particular, we do not need any assumed bound on the solution to 
get the error estimate. 

The mathematical analysis of numerical schemes for the discretization of the steady and/or non steady 
compressible Navier-Stokes and/or compressible Stokes equations has been the object of some recent 
works. The convergence of the discrete solutions to the weak solutions of the compressible stationary 
Stokes was shown for a finite volume- non conforming PI finite element mannnB] and for the well- 
known MAC scheme which was introduced in m and is widely used in computational fluid dynamics 
(see e.g. [29]). The unsteady Stokes problem was also discretized by some other discretization schemes on 
a reformulation of the problem, which were proven to be convergent [241426] . The unsteady barotropic 
Navier-Stokes equations was recently investigated in m in the case 7 > 3 (there is a real difficulty in 
the realistic case 7 < 3 arising from the treatment of the non linear convective term). However, in these 
works, the rate of convergence is not provided; in fact, to the best of our knowledge, no error analysis 
has yet been performed for any of the numerical schemes that have been designed for the compressible 
Navier-Stokes equations, in spite of its great importance for the numerical analysis of the equations and 
for the mathematical simulations of compressible fluid flows. We present here a general technique to 
obtain an error analysis and apply it to one of the available numerical schemes. To the best of our 
knowledge, this is the first result of this type in the mathematical literature on the subject. 

To achieve the goal, we systematically use the relative energy method on the discrete level. From 
this point of view, this paper is as valuable for the introduced methodology as for the result itself. Here, 
we apply the method to the scheme of m- In spite of the fact that this latter scheme is not used in 
practice (see e.g. [28] for a related schemes used in industrial codes), we begin the error analysis with 
the scheme m because of its readily available convergence proof. In fact, we aim to use this approach 
to investigate the numerical errors of less academic numerical schemes, such as the finite volume - non 
conforming PI finite element [171121711211123] or the MAC scheme mm- 

The paper is organized as follows. After recalling the fundamental setting of the problem and the 
relative energy inequality in the continuous case in Section [2] we proceed in Section [3] to the discretiza¬ 
tion: we introduce the discrete functional spaces and the definition of the numerical scheme, and state 
the main result of the paper, that is the error estimate formulated in Theorem 13.11 The remaining 
sections are devoted to the proof of Theorem 13.11 

• In Section [4] we recall the existence theorem for the numerical scheme (Lemma 13.11) and derive 
estimates provided by the scheme. 

• In Sectional we derive the discrete intrinsic version of the relative energy inequality for the solutions 
of the numerical scheme (see Theorem ED. 
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• The relative energy inequality is transformed to a more convenient form in Section [ 6 j see Lemma 

P 

• Finally, in Section [71 we investigate the form of the discrete relative energy inequality with the 
test function being a strong solution to the original problem. This investigation is formulated in 
Lemma Q and finally leads to a Gronwall type estimate formulated in Lemma 18.11 The latter 
yields the error estimates and finishes the proof of the main result. 

Fundamental properties of the discrete functional spaces needed throughout the paper are reported in 
Appendix (Section ED . Some of them (especially those referring to the L p setting, p / 2 that are not 
currently available in the mathematical literature) are proved. Section E] is therefore of the independent 
interest. 


2 - The continuous problem 


The aim of this section is to recall some fundamental notions and results. We begin by the definition of 
weak solutions to problem (ED- C3D- 

Definition 2.1 (Weak solutions). Let go '■ Ll ^ [0,+oo) and uq : 0 —>• W l with finite energy Eq = 
/n(l£o| u o | 2 + H(qq)) dx and finite mass 0 < Ma = f Q go dx. We shall say that the pair (g, u ) is a weak 
solution to the problem CCD (USD emanating from the initial data (go,tio) if: 

(a) g € L°°(0, T; L 1 (f2)), g> 0 a.e. in (0 ,T), undue L 2 (0, T; W 0 1,2 (P)). 

(b) g e C we ak([0, T]; L 1 (0)), and the continuity equation (II.lab is satisfied in the following weak sense 

J gipdx =J J (gd t p + gu-Vx<p) dxdt, Vr G [0,T], Vy? € C c °°([0, T] x 0). (2.1) 

(c) gu e C we ak([0, T] ; L 1 (Q)), and the momentum equation (11.161) is satisfied in the weak sense, 

J £m-<p?dx| = J J (^gu ■ dt^p + gu ® u : V</? + p(g) div<^ dx dt 

— J J ^Vm : \7 x (pdxdt+(ti + A)divudiv<^j dxdt, Vr € [0,T], Vip e C“([0,T] x 0;M 3 ). 

( 2 . 2 ) 


(d) The following energy inequality is satisfied 


[ (~g\u\ 2 + H(g)\ dx + [ [ (fi\ Vn | 2 + (fj, + A)| divn| 2 ) dx dt < 0, for a.a. tG(0,T), (2.3) 
Jn v 2 > o Jo Jn v ' 


with H(g) = 



(2.4) 


Here and hereafter the symbol / g dx is meant for / g(r, x) dx — / go{x)dx. 

Jn Jn Jn 

In the above definition, we tacitly assume that all the integrals in the formulas CEBH23D are defined 
and we recall that C wea k([0, T]; L 1 (0)) is the space of functions of L°°([0, T]\ L l (fil)) which are continuous 
for the weak topology. 

We notice that the function g i-> H(g) is a solution of the ordinary differential equation gH'(g) — 
H(g) = p(g) with the constant of integration fixed such that H( 1) = 0. 
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Note that the existence of weak solutions emanating from the finite energy initial data is well-known 
on bounded Lipschitz domains under assumptions (11.51) and (11.61) provided 7 > d/(d — 1), see Lions (30J 
for "large” values of 7 , Feireisl and coauthors [H] for 7 > d/(d — 1). 

Let us now introduce the notion of relative energy. We first introduce the function 


E : [0, 00 ) X (0, 00 ) —> M, 

(g, r) 1-7 E(g\r) = H(g) — H'(r)(g — r) — H(r ), 


(2.5) 


where H is defined by (12.41) . Due to the monotonicity hypothesis in (11.51) . H is strictly convex on [0, 00 ), 
and therefore 

E(g\r) > 0 and E(g\r ) =0 47 g = r. 

In order to measure a “distance” between a weak solution (g, u ) of the compressible Navier-Stokes system 
and any other state (r, U) of the fluid , we introduce the relative energy functional, defined by 


£(g, u r,U) = j {^g\ u ~ U\ 2 + E(g | r)j dx. 


( 2 . 6 ) 


It was proved recently in m that, provided assumption (ESI) holds, any weak solution satisfies the 
following so-called relative energy inequality 


£ 


(g,u r,U^ (t) — £ (g,u r,U^j (0) + j J (/r|V(w - U)\ 2 + (g + A)| div(u - £/)| 2 ) dxdt 

~ Jo J (/ iVi7: v (^-w) + (/r + A) div D div (I7-w)) dx dt 

+ / [ gdtU ■ (U — u) dxdt+ [ [ gu-'S/U ■ (U — u) dxdt 

Jo Jn Jo J n 

p(g)divU dxdt+ / / (r — g)dtH'{r) dxdt — / / gVH'(r)-u dxdt 

Jo Jn Jo Jn 


(2.7) 


/0 Jn 


for a.a. r G (0, T), and for any pair of test functions 

r G C' 1 ([0,T] xfi), r > 0, UgC\[0,T] xD;M 3 ), U | an = 0. 


The stability of strong solutions in the class of weak solutions is stated in the following proposition. 

Proposition 2.1 (Estimate on the relative energy). Let Q be a Lipschitz domain. Assume that the 
viscosity coefficients satisfy assumptions {Z3])> that the pressure p is a twice continuously differentiable 
function on (0, 00 ) satisfying \1.5\) and <11. 61) . and that (g, u) is a weak solution to problem (11.11) - (11.31) 
emanating from initial data (go > 0 ,uq), with finite energy Eq and finite mass Mo = godx > 0. Let 
(r,U) in the class 

! r G C 1 ([0, T] x fl), 0 < r = min_ r(t, x) < r(t, x) <r = max r(t,x), 

{t,X)(fzQrp (t)X)£.Qrp 

( 2 . 8 ) 

[/GC^O.TlxlI;! 3 ), t/|( 0 ,T)xSn = 0 

be a (strong) solution of problem (11,11) emanating from the initial data (ro,Uo). Then there exists 
c = c(T,n,M 0 ,Eo,r,r, \p'\ c i(]r,r]), ||(Vr, d t r, U, VC/, 9 t C7)|| L oo ( Q T;R i9 ) ) > 0 
such that for almost all t G (0,T), 


£(g,u r, U)(t) < c£(g 0 ,u 0 r 0 ,U 0 ). 


(2.9) 


This estimate (implying among others the weak-strong uniqueness) was proved in [L2] (see also [T5]) 
for pressure laws Ea with 7 > d/(d — 1). It remains valid under weaker hypothesis on the pressure, 
such as (11.61) with 7 > 1 ; this can be proved using ideas introduced in | 2 ] and m- 
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3 - The numerical scheme 


3.1 Partition of the domain 

We suppose that P is a bounded domain of M d , polygonal if d = 2 and polyhedral if d = 3. Let T be a 
decomposition of the domain P in tetrahedra, which we call hereafter a triangulation of P, regardless of 
the space dimension. By S(K), we denote the set of the edges (d = 2) or faces (d = 3) a of the element 
AT G T called hereafter faces, regardless of the dimension. The set of all faces of the mesh is denoted by 
£; the set of faces included in the boundary <9P of P is denoted by £ ex t and the set of internal faces (i.e 
£ \ £ ex t) is denoted by £; n t• The triangulation T is assumed to be regular in the usual sense of the finite 
element literature (see e.g. 0 ), and in particular, T satisfies the following properties: 

• P = 

• if (. K,L ) € T 2 , then KOL = 0 or K n L is a vertex or K n L is a common face of K and A; in the 
latter case it is denoted by A"|L. 

For each internal face of the mesh a = K\L, stands for the normal vector of a, oriented from 

K to L (so that n„j<c = — n CT ./-J. We denote by \K\ and |cr| the (d and d — 1 dimensional) Lebesgue 
measure of the tetrahedron AT and of the face o respectively, and by hx and h a the diameter of K and 
a respectively. We measure the regularity of the mesh thanks to the parameter 6 defined by 

6 = inf{|—, K G T} (3.1) 

h K 

where £k stands for the diameter of the largest ball included in 
the maximal size of the mesh, 

h = max his. 

KeT 

The triangulation T is said to be regular if it satisfies 

6>6 o > 0. (3.3) 


K. Last but not least we denote by h 

(3.2) 


3.2 Discrete function spaces 

Let T be a mesh of P. We denote by L^(P) the space of piecewise constant functions on the cells of 
the mesh; the space L) l (P) is the approximation space for the pressure and density. For 1 < p < oo, the 
mapping 

«« MiJ(n) = II?IIl, (1! ) = ( E \K\\<iK\ r ) llV 

KeT 

is a norm on L/j(P). We also introduce spaces of non-negative and positive functions: 

L+(P) = {q € L h { P), q K > 0, VW G T}, L++(P) = {q € L h { P), q K > 0, VK € T}. 

The approximation space for the velocity field is the space W/ l (P) = 14(P;M d ), where 14,(P) is the non 
conforming piecewise linear finite element space 00 defined by:. 

Vfc(P) = {u € L 2 (P), VK G r, V\ K G Pi (AT), 

Vex G £\nt , o' = AT|L, f v\ K dS = f v\ L dS, Mo G £ ext , [ vdS = 0}, (3.4) 

J (7 J (7 J a 

where Pi (A") denotes the space of affine functions on K and dS the integration with respect to the 
(d — l)-dimensional Lebesgue measure on the face a. Each element v G V4(P) can be written in the 
form 

v(x) = ^2 x € P ( 3 - 5 ) 
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where the set {(Po}o££ int C I4(fi) is the classical basis determined by 


V(ct,ct') E £? nt , f 

J a 


7>C 


dS = 5„ 


Va' € S. 


ext 5 


Vt 


d5 = 0 


( 3 . 6 ) 


and {fo-jo-efint C R is the set of degrees of freedom relative to v. Notice that 14(0) approximates the 
functions with zero traces in the sense that for all elements in 14(0), v a = 0 provided a E £ ex t- Since 
only the continuity of the integral over each face of the mesh is imposed, the functions in 14 ( 0 ) may be 
discontinuous through each face; the discretization is thus nonconforming in W 1,p (0;M d ), 1 < p < oo. 
Finally, we notice that for any 1 < p < oo the expression 


V \v£(n) ( 

KeT 


l/p 


is a norm on 14(0) and we denote by VF (O) the space 14(0) endowed with this norm. 

We finish this section by introducing some notations. For a function v in L 1 (0), we set 



( 3 . 7 ) 


so that v E L/j(O). Here and in what follows, Ik is the characteristic function of K. 
If v € bF 1,p (0), we set 

v a = -—- / vdS for er E £. 

\a\ Jo 


Finally, if v E Wo’ P (0), we set 


Vh(x) = v a ip a (x), x E 0. 


( 3 . 8 ) 

( 3 . 9 ) 


so that Vh £ 14(0). In accordance with the above notation, for v E Wq’ p { O), the symbol Vh means 
Vh{x) = E CT e£int Vo4>o(x), and the symbol v h , K = yk] Ik v h(x)dx and the symbol v^ a = [(u/4]" p . 


3.3 Discrete equations 

Let us consider a partition 0 = to < ti < ••• < tjv = T of the time interval [0, T], which, for the sake 
of simplicity, we suppose uniform. Let k be the constant time step k = t n — t n — i for n = 1, The 

density field g(t n ,x ) and the velocity field u(t n ,x) will be approximated by the quantities 

e n (x) = QkIk(x), u n (x) = £ Kvo(x), ( 3 . 10 ) 

KeT oes 

where the approximate densities (^)iteT,n=i,...,iv and velocities (u™) a e£ int ,n=i,...,JV are the discrete un¬ 
knowns (with E M + and ti” E R rf ). 

For the future convenience, we denote here and hereafter, 


N 


N 


8(t,x) = 0 n (*) 1 [n-l,n)W» “(*,®) = E , * n ( iI ) 1 [n-l.n)( t ) 


( 3 . 11 ) 


n =1 n=l 

and recall that the usual Lebesgue norms of these functions read 


|L°°(0,T;LP(n) = n ™ax^ ||0 ||z,P(U), ll' u lll,P(0,T;L9(n;R 3 ) 


N 

= *(£ 

n=l 


1 iiL9(n ; R 3 ) 


1/p 


( 3 . 12 ) 


Starting from this point, unlike in Section [I] here and hereafter, the couple (g, u) respectively (g n , u n ) 
introduced in (I3.10I1XTID denote always exclusively a discrete numerical solution. 
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The numerical scheme consists in writing the equations that are solved to determine these discrete 
unknowns. In order to ensure the positivity of the approximate densities, we shall use an upwinding 
technique for the density in the mass equation. For q E L/^fl) and u E W^fi), the upwinding of q with 
respect to u is defined, for a = K\L E £; n t by 


q 


up _ 
a 


q K if u a ■ n„, K > 0 
q L if u a ■ n^ K < 0 , 


(3.13) 


so that 

Y qa P u cr ■ n a>K = Y ( QkIuv ■ n CT) ^] + - q L [u a ■ n ^]"), 

ae£(K) ae£(K) 

rj=K\L 

where a + = max(a,0), a~ = — min(a, 0). 

Let us then consider the following numerical scheme m 

Given ( g°,u °) € L£(fl) x W h (p) find (g n )i< n <N C {L h {Q)) N , (u n ) i< n <jv C such that for 

all n = 1,..., N 


n n _ n n- 1 

\K\ Qk Q * + E M#" p 

*e£{K) 

[K ■ rifj K ] = 0 , 

VI< E T, 

(3.14a) 

E - e* 1 ”* 1 ) ■ 

K&T H 

v k+ Y 

\a\gY P K’ up i< ■ ™*,k] ■ v K 


KeT aES(K) 




- E ?«> Yj H Vg- ■ tier.ft' + /i ^ / Vit n : Vu dx (3.14b) 

A'eT o-g£(a:) k&t Jk 

+ (fi + \) Y ( divi/Miv-u dx = 0, Vi> E 
K&T K 

Note that the boundary condition tt” = 0 if cr E £ ex t is ensured by the definition of the space Vf-fiii). 
Note also that if a E £T int , cr = K\L, one has, following (13.71) and (13.131) . 

u"’ up _ u n _ / u n^ x j ( ^ x j£ u'f . n a K >0 and ■u” ,up = uf = t~t [ u n {x) dx if w" • n a k < 0 . 

\K\ Jk ’ \L\ Jl 

It is well known that any solution ( Q n )i< n <N C ( Lh(Q)) N satisfies g n > 0 thanks to the upwind 
choice in (13.14aD (see e.g. [171127]). Furthermore, summing (I3.14al) over K E T immediately yields the 
total conservation of mass, which reads: 

Vn = 1, ...N, [ g n dx = [ g° dx. (3.15) 

Jfl J£l 

We finally state in this section the existence result, which can be proved by a topological degree 
argument, mm- 

Proposition 3.1 (Existence). Let (g°, n°)E L)J) + (I1) x Whip). Under assumptions (11.41) and (USD, 
Problem (|3.14|) admits at least one solution 

(g n )l<n<N € [Lt + m N ,(.V n )l<n<N € [W^fi)]*. 
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3.4 Main result: error estimate 


Let (r,U) :[0,T]x(!h> (0, oo) x R 3 be C 2 functions such that U = 0 on d£l. Let ( g,u ) be a solution 
of the discrete problem (13.1411 . Inspired by (12.fill , we introduce the discrete relative energy functional 


£(f? n X 


r n ,U n ) = 


KG' 


g n \u n -Uff+ E{g n \r n ) dx 


(3.16) 


y. 1^ ~ ^h,A'| 2 + ^(^A'l r A'))i 


AeT 


where 


r n (x ) = r(t n ,a;), £/ n (x) = U(t n ,x), n = 0,...,N, 

( g n ,u n ) is defined in (13.1011 . and if is defined by (12.511 . Let us finally introduce the notations 

M 0 = y \K\g%, and Eq = £ |tf| {\q°k\u ° k \ 2 + #(&)). 


(3.17) 


Ae/c 


A e/C 


Now, we are ready to state the main result of this paper. For the sake of clarity, we shall state the 
theorem and perform the proofs only in the most interesting three dimensional case. The modifications 
to be done for the two dimensional case, which is in fact more simple, are mostly due to the different 
Sobolev embedings, and are left to the interested reader. 

Theorem 3.1 (Error estimate). Let 9 q > 0 and T be a regular triangulation of a bounded polyhedral 
domain £l C R 3 introduced in Section no such that 9 > 9 q, where 9 is defined in (EU). Let p be a twice 
continuously differentiable function satisfying assumptions with 7 > 3/2, and the additional 

assumption m in the case 7 < 2. Let the viscosity coefficients satisfy assumptions (Q- Suppose that 
(g°,u°) € Lff(Ll) x Wh(Ll) and that ( g n )i< n <N C [L^(Q)] N , (u n )i< n <jv C [W^fi)] - ^ is a solution of 
the discrete problem (13.1411 . Let (r,U) in the class 


r G C 2 ([0, T] x fl), 0 < r := min_ <r(t,x)<r:= max r(t,x), 




(t,x)£Q T 


(3.18a) 

(3.18b) 


U € C 2 ([0,T] xll;R 3 ), U\ dn = 0 

be a (strong) solution of problem (11.111 . Then there exists 

c = c(t, |n|,diam(0),6»o,7,Mo,Eo,r,r, 

bVfcr]), \\(Vr,d t r,d t Vr,d 2 r,UyU,W 2 U,d t U,d 2 UAyU)\\ Loo{QT . R6S) ) € (0,+ 00 ) 
(independent of h, k) such that for any m = 1,... ,N, 

m ,, 

£(g m ,u m \r m ,U m )+kJ2J2 / |V x (w n - U/)| 2 dx < c{s(g°, u°|r°, U°) + h A + Vk) , (3.19) 

n=l AeT ^ K 


where 

A= [^ 1/7 € (3/2,2], 

11/2 ifl>2. 


(3.20) 


Starting from this point, unlike in Section [T] here and hereafter, the symbol £ refers always to the 
discrete relative energy functional defined in (13.1611 . 
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Remark 3.1. 

Assumptions \d.l8 1) on the regularity of the strong solution (r, U) in Theorem ld. 1\ may he slightly relaxed: 
It is enough to suppose 

(r, U) € C' 1 ([0, T] x O; M 4 ), V 2 U € C([0, T] x H; R 3 ), 0 < inf_ r(t,x), 

( t,x)eQ T 


d$r € L 1 (0,T-,L'i'(SI)), d t Vr £ L 2 (0,T-,L^/^~ 6 \n-,R 3 )), (d 2 U, 3 t VU) € L 2 (0, T; L 6 / 5 (fl; R 12 )). 

The constant in the error estimate depends on r and the norms of r and U in these spaces. This 
improvement is at the price of more technicalities in estimates of several residual terms, namely in 
estimates \6.d\ - \6. 5lJ . {6. Iff) , \ 6. 22\) . 17. .hi) , i l 7. Ill F T . ill) and &8.3I) . 

Remark 3.2. 


1. Theorem VI. 11 holds also for two dimensional bounded polygonal domains under the assumption that 
7 > 1. Assumption on the asymptotic behavior of pressure near 0 is no more necessary in 

this case. The value of A in the error estimate id. 191) is 


A = 


^ i/7 € (1,2], 

1 if 7 > 2 . 


2. Suppose that the discrete initial data (g°, u°) coincide with the projection (f°,U°) of the initial 
data determining the strong solution. Then formula id. 191) provides in terms of classical Lebesgue 
spaces the following bounds: 


|| m _ m II2 

lit' ' llL 2 (Qn{r/2< e m <2r}) 


^ U \\L' 2 {n^{r/2<e rn <2r}) — C (j l + 


for the "essential part" of the solution (where the numerical density remains bounded from above 
and from below outside zero), and 

\{s m < L/2)\ + \{e m > 2r}| + lli>”III, (!! n( c ~> 2 r}) + 11^1^” " U"f || 11(nn ( c ».> 2 r)) < c(h A + Vk) 

for the "residual part" of the solution, where the numerical density can be "close" to zero or infinity. 
(In the above formula, for B C Tl, \B\ denotes the Lebesgue measure of B.) 

Moreover, in the particular case of p(g) = g 2 (that however represents a non physical situation) 
E(g\r) = (g — r) 2 and the error estimate id. 191) reads 


Q m ~ r m \\h m + ll£ m |u m - U"f || Ll(a) < c(Vh + Vk) 


3. Theorem ld.ll can be viewed as a discrete version of Proposition 12.11 It is to be noticed that the 
assumptions on the constitutive law for pressure guaranteeing the error estimates for the scheme 
(13.141) are somewhat stronger (7 > 3/2j than the assumptions needed for the stability in the contin¬ 
uous case ( r y>l). The threshold value 7 = 3/2 is however in accordance with the existence theory 
of weak solutions. The assumptions on the regularity of the strong solution to be compared with 
the discrete solution in the scheme are slightly stronger than those needed to establish the stability 
estimates in the continuous case. 

4- If d = 3, we notice that the assumptions on the pressure (as function of the density) in Theorem 
F P are compatible with the isentropic case p(g) = g 7 for all values 7 > 3/2. 


10 



















5. The scheme JEHU contains in addition artificial stabilizing terms both in the continuity and momen¬ 
tum equations. These terms are necessary for the convergence proof in even for the large values 
of 7 . It is to be noticed that the error estimate in Theorem \S.1\ is formulated for the numerical 
scheme without these stabilizing terms. Of course similar error estimate is a fortiori valid also for 
the scheme with the stabilizing terms, however, this issue is not discussed in the present paper. 

The rest of the paper is devoted to the proof of Theorem 13.11 For the sake of simplicity, and in 
order to simplify notation, we present the proof for the uniformly regular mesh meaning that there exist 
positive numbers q = Cj(#o) such that 

C] Iik <h< c 2 h a < c 3 h,K, Cl \K\ < \cr\h < c 2 \a\hK < c 3 \a\h a < c±\K\ (3.21) 

for any K E T and any a E £. The necessary (small) modifications needed to accommodate the regular 
mesh satisfying only (13.31) are straightforward. Even with this simplification the proof is quite involved, 
and some details have to be necessarily omitted to keep its length within reasonable bounds. The reader 
can eventually find them in the extended version of this paper available on ArXiv m- 


4 - Mesh independent estimates 

We start by a remark on the notation. From now on, the letter c denotes positive numbers that may 
tacitly depend on T, |fi|, diarn(Q), 7 , a, do, A and /r, and on other parameters; the dependency on 
these other parameters (if any) is always explicitly indicated in the arguments of these numbers. These 
numbers can take different values even in the same formula. They are always independent of the size of 
the discretisation k and h. 


4.1 Energy Identity 

Our analysis starts with an energy inequality, which is crucial both in the convergence analysis and in 
the error analysis. We recall this energy estimate which is already given in EZ3, along with its proof for 
the sake of completeness. 

Lemma 4.1. Let ( g°,u° ) E L~ } J)(P) x Wf,(Q) and suppose that ( g n )\< n <N £ )l<n<AT £ 

mm N is a solution of the discrete problem ([3.141) with the pressure p satisfying condition 111.5)) . Then 
there exist 

Qa £ [min(^, gf), max(^, gf)}, a = K\L E £ int , n = l,...,N 
—n—i, n € ^^ max (^- 1 ) 0 ^)], K E T, n = 1,..., N 

such that 


y . \^\( < 2^k\ u 'k\ 2 + ^{q'k)) ~ ^2 \I x \{f58%\ u %V + H {q°k) 


KeT 


K&T 


m ,, 

+ k (jjl / iV^n ”! 2 dx + (g + A) / |divu n | 2 dx^ 

n =l KeT V ^ K ^ K 

+ Wif 1 ] + [D t Ti atl ] + re£“'] + = 0, (4.1) 
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for all m = 1,..., N, where 


[A?i A “'l = E E \K\el - 1 


.n —1 |2 


n=i /feT 

771 


[A 7 n!e Ael ] = E E 


I^a: 8k 


71 — 1 I 2 


-K 


7i=i KeT 

771 


[oS A “'] = tE E 


I I n un ( U K ~ U L? I n I 

klgg’ p 2 K ' n ffl A-|, 


n=lcr=A'|LGfint 

771 


[o,”^'] = *E E M ff "® 

”=1 a=K\L££ in t 


i--- • n CT A-1. 


(4.2a) 


(4.2b) 


(4.2c) 


(4.2d) 


Proof. Mimicking the formal derivation of the total energy conservation in the continuous case we take 
as test function v = u n in the discrete momentum equation (|3.14bl) n and obtain 


I\ + I2 + -^3 + ^4 — 0 , 


(4.3) 


where 

h = E ¥<« - eT'u"-'' 


A'eT 


k 


,< j ■ «&, '2 = E E kk: “ p «’ up ■ <■ «• 

KeT aeS(K) 
a=K\L 


h = ~ E E |cb(gA:)KJ ■ n <r,A:], h = E / : Vi / 1 + (// + A)divu n divu") dx. 

KeTae£(K ) ^^ J K 


KeT 


t=K\L 


Next, we multiply the continuity equation (I3.14all j r by and sum over all K € T. We get 

^5 + ^6 = 0 

with / 5 = - £ - &-‘)l<-| 2 and /a = - E E 5M<?r p [< ' >V*]I“*I 2 - 


( 4 . 4 ) 


ATeT 


A'eT ae£(K) 
a=K\L 


Finally, we multiply the continuity equation (13.14aD W by H’(o n K ) and sum over all K £ T. We obtain 


h + h — 0, 

with Ir = Y, ^(Qk - sTTm and J 8 = E E M# Up K • n a , K ]H\e n K ). 

KeT KeToe£(K) 

a=K\L 

We now sum formulas (14.311 — (14.511 in several steps. 

Step 1 : Term I\ + I 7 . We verify by a direct calculation that 


(4.5) 


h = e ^(5 a wi 2 - hrvrT) + e t 1 *' 

KeT L 1 1 KeT k 

In order to transform the term It, we employ the Taylor formula 

H'm(& - A" 1 ) = urn - H{elT + iff"@r‘’”)« - A -1 ) 2 , 


1 


I A 1 n-l\ U K - W A rl | 2 
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where g^ 1,n £ [min (g 7 ^ 1 ,^-),max(^ 1 , 6k)\- Consequently, 


n— 1 


n— 1 


h + h=Y, ^(5*K-l 2 - Lr‘KT 2 ) + £ x(»(<£) - T*" 1 )) 


Tl 


AeT 


AeT 

+ e Z ** 1 ' 2 + £ 7 * 1 ' 2 - ( 4 . 6 ) 


k 


KeT KeT 

Step 2 : Term I 2 + Iq- The contribution of the face o = K\L to the sum I 2 + Iq reads, by virtue of 


0 ), 


M [< • n <x,A] + £a(W| 2 - U K ■ U L - \\ U k \ 2 + tKI 2 ) 


+ |(t| k • ^,l]+ ql (l^L [ 2 - < ■ < - ^Kl 2 + ^Kl 2 


Consequently, 

^2 + / 6 = E kl|<-n^| g r p (< ~^ )2 - 

0~ = K | -£/££int 

Step 3 : Term T 3 + /§ • We have 


(4.7) 


*=E E 


cr • n, 


AeT ctg£(K) 
a=K\L 


,a] (i7'(^)(^ up -^) + 77(^)) 


+ E E M K • 

Kdz'T a(E:£(K) 
a=K\L 


,k] (<£»'(<&) - H(e%) 


Recalling (13.131) . we may write the contribution of the face cr = K\L to the first sum in Is; it reads 

M K • ^, K ]+ (H(e n K ) - n'{el)(e n K - el) - 

+ M K ■ >vzJ + (1H(el) - HXeDlel - el) - ff(efc)). 

Recalling that rH'(r) — H(r) = p(r), we get, employing the Taylor formula 

\ 2 

(g n _ Q n 

h + h = E \<-n a , K \H"(gZ)- K L ‘ 
o-=A|Le£mt 

with some ~g£ € [min(^, g"),max(g|, £>2)]. 

Step 4: Conclusion 

Collecting the results of Steps 1-3 we arrive at 


£vf( 


AeT 


,Qk\ u k\ 2 - 6k 1 \ u k E) + E H (8k *)) + E (t ^ \V x u n \ 2 dx 


AeT 


AeT 


+ o*+a)/ |div U -i 2 dd + £ l|W l W ' E' 2 + £ 

4a ' TtL k 2 k 2 


„n _n—112 


AeT 


AeT 


cre£int 

<T=A|L 


1 I I n,up ( U K ~~ u l) Z i n 1 , |_|z j" ( -n 

+ 2^ l°W- 2 -I™<t- j Va|+ 2^ l 0 ^ 


(, Qk-Ql 


\ u a ‘ n cr,A| = 0. (4.8) 


C(z£int 

a=K\L 


At this stage, we get the statement of Lemma 14.11 by multiplying (I4.8l) n by k and summing from n = 1 
to n = m. Lemma 14.11 is proved. □ 
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4.2 Estimates 

We have the following corollary of Lemma 14.11 


Corollary 4.1. (1) Under assumptions of Lemma there exists c 

h and k) such that 

l W lz, 2 (0,T;V^(n;R 3 ) — C 
IMlL 2 (0,T;L 6 (n ; M 3 )) — C 

II /v 2 ii 

\\Q U IIl oo ( 0,T;L 1 (0)) < c - 

(2) If in addition the pressure satisfies assumption H1.6\) then 


c(Mq,Eo) > 0 (independent of 

(4.9) 

(4.10) 

(4.11) 


\\q\\l^(o,t-,l^q)) < c 


(4.12) 


(3) If the pair (r, U) belongs to the class \3.18\) there exists c = c(Mq, Eo,r , r, ||C7, VL/|| LO o(q t . r i 2 )) > 0 
such that for all n = 1,... , N, 


£(g n ,u n \r n ,U n ) < c, 

where the discrete relative energy £ is defined in \3.1(i\) . 


(4.13) 


Proof. Recall that 


' U Il 2 (0,T;V 2 (£‘2;R 3 ) 



X7 x u n \ 2 dx; 


the estimate (14.91) follows from (14.11) . The estimate (14.101) holds due to imbedding (19.301) in Lemma 19.31 
and bound (14.91) . The estimate (14.111) is just a short transcription of the bound for the kinetic energy in 

(ED- 

We prove estimate (14.121) . First, we deduce from (11.51) and the definition (12.41) of H that 0 < —H(g) < 
c\ with some c\ > 0, provided 0 < g < 1 and H{g) > 0 if g > 1. This fact in combination with the 
bound for f n H(g)dx derived in (14.11) yields 


/ \H(g)\dx < c < oo. 
Jn 


(4.14) 


Second, relations (11.5111.71) imply that there are g > 1 and 0 < p < p < oo such that 


f Q a Po/2 < ^if o < 0 < 1/e, 1 
l p < < p if i/e < e < e, \ 

{ ^ 7_2 Poo/2 < ^ if 6 > 2g J 

Using these bounds and the definition (El of H we verify that 


Q 7 < c(\H(g)\ + g+ 1) 


with a convenient positive constant c. Now, bound (14.121) follows readily from the boundedness of 
f n g m dx = Jf,KeT\K \&k an< i fnH{Q m )d x = J2xeT\K\H(Q k) established in (13.151) and (14.11) . 

Finally, to get (14.131) . we have employed (12.51) . (I3.1(il) . (13.151) . (14.141) to estimate f Q E{g n \f n ) dx and 
(|4.11|) . (19.31) . (I9.22|) to evaluate J2k&t Ik 0K\^h,K ~ u k\ 2 d x■ □ 

The following estimates are obtained thanks to the numerical diffusion due to the upwinding, as is 
classical in the framework of hyperbolic conservation laws, see e.g. [9J. 
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Lemma 4.2 (Dissipation estimates on the density). Let ( g°,u °) € L^(fi) x Wh{Ll). Suppose that 
(Q n )i<n<N C (■u n )i< n <iv C [W) l ] Ar (D) is a solution of problem (I3.14|) . Finally assume that 

the pressure satisfies hypotheses (ESI) and (EJh Then we have: 

(1) If 7 > 2 then there exists c = 0 ( 7 , 9q,Eq) > 0 such that 


N 


fc E E 

n ~ 1 cr—|LE£int 


(el ~ el? 

max(g^, off) 


uf ■ n^K'l < c. 


(4.15) 


(2) 1/7 € [1, 2) and the pressure satisfies additionally assumption hi. 71) then there exists c = c(Mo, 77o) > 
0 snch that 

N ( n n — n n 'l 2 

fc E E l 0- !?—— «!i2- 7 1 {gg>i} l u " • n ^l 


T- 1 (J— K\LG£int 


N 


+ k E E M(eSr - ^l) 2 i {?S<i} K • n <x,/d < c, (4.16) 

ti= 1 <j=K\L€.£i n t 

where the numbers ~q% are defined in Lemma o 

Proof. We start by proving the simpler statement (2). Taking into account the continuity of the pressure, 
we deduce from assumptions (11.61) and (11.71) that there exist numbers p 0 > 0 , > 0 such that 


H"(s) > 


[jjSy, ifs>l, 

[p 0 s a > p Q , if s < 1 , 


whence, splitting the sum in the definition of the term (see (14.2dl) f into two sums, where (a, n) 

satisfies > 1 for the first one and ~gf < 1 for the second, we obtain the desired result. 

Let us now turn to the proof of statement (1). Multiplying the discrete continuity equation (13. 14a 1) 7- 
by In q\ and summing over K £ T, we get 

r.n _ n n ~ 1 

E \ R \ - k K '°a+ E E (taes-)sr p <-^,jf = o. 


AeT 


K&T a££(K),cr=K\L 

By virtue of the convexity of the function g 1 —?• £> In g — g on the positive real line, and due to the Taylor 
formula, we have 

Q H K In Ql ~ (Tk 1 ^ Q U k X ~ (el - &') < In el(el ~ Q^)\ 
whence, thanks to the mass conservation (13.151) and the definition of g'fP, we arrive at 

Ji— 1 1 - n—1 

k 


E \K\ IkhIk esr’i^ggr 1 + y k |* K . n „ iA . l+ (i n *_i n£ , 2 ) 


AeT 


0"€:£int 

a=K\L 


+ E l c7 I^L[<-^,L] + (ln^- 1 n^) <0, 


O'Gfint 
cr=A|A 


or equivalently 

7 £ |a|K-n^]+(^(ln^-ln^)-(^-^))+fc £ \a\[u^ L } + (gl(ln gf-ln g n K )-(g n L -g n K ) 


cGfint 
ct=A|A 


O'S^int 

<r=A|L 


AeT 


O'eTnt 
(7=A | L 


- £ \K\ (g n K In el - g^- 1 In gffi 1 ) + k £ M ([< ■ ^,a] + (^ - &)) + K ■ n a , L ] + (gl - g n L ))] 

(4.17) 
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From m Lemma C.5], we know that if ip and if are functions in C 1 ((0, oo); R) such that sif'(s) = ip'(s) 
for all s E (0, oo), then for any (a, b ) E (0, oo) 2 there exits c € [a, b] such that 


1 


0(6) - if(a))b - (ip(b) - <p(a)) = f-(b - a) V(c). 


Applying this result with if(s) = In s, (p(s) = s we obtain that the left hand side of (14.1711 is greater or 
equal to 


0"€=£int 

a=K\L 


On the other hand, the first term at the right hand side is bounded from above by Finally 

the second term at the right hand side is equal to 

~ k Y L 8K di ™ n < k ]T II^A'||L2(A-)l|divu n || Z/ 2 (K) , 


A'eT 


IK 


K&T 


whence bounded from above by k\\u n \\ V 2^. R 3' ) \\g n \\L2^, where we have used the Holder inequality and 
the definition of the V) 2 (Q)-norm. The statement (1) of Lemma 14.21 now follows from the estimates of 
Corollary 14.11 

□ 


5 - Exact relative energy inequality for the discrete problem 


The goal of this section is to prove the discrete version of the relative energy inequality. 

Theorem 5.1. Suppose that O C M 3 is a polyhedral domain and T its regular triangulation introduced 
in Section UTH Let p satisfy hypotheses J 1. 5\) and the viscosity coefficient p, A obey P^\ >- Let (g°, u°) E 
L^(fi) x Wh(fl) and suppose that ( g n )i< n <N € [£^(17)]^, ( u n )i< n <N € [TV/ l (r2)] iV is a solution of the 
discrete problem (13.141) . Then there holds for all m = 1,..., N, 


K&T 


Y. - U™ K \ 2 - Qk\ u K - Uy l K \ 2 ') + Y, l-fr I ^(Qk^k) — ^(^A'l r A:)) 

AeT 

772 r. r. 

+ k Y Y (/i \S7 x {u n - Uf) | 2 dx + (fi + A) |di v(u n - Uf)\ 2 dx) 

n=lK&T JK JK 

m „ r. 

E E (W : V*(E# - “") dx + (g + X) divl/?div(E# - u n ) dx) 

ntlAET JK JK ’ 

m TT n TT n ~ 1 TT n ~^ i rrn 

+k ee — k ~ - u * 1 


,U? k + ULt s 

.72,up ( n,ij ^.n,up\ 1 - 7-71 \„.n _ , 

’ { -2- a ) ' U h,K[ u <y • n ^A-J 


ra=l AeT 

m 

-^e e y \ a \e< 

72=1 KeT cre£(K) 
a=K\L 

m 

-‘EE E W\p(eMuK-n„, K ] 

n =1 KeT ae£(K) 
cr=K\L 

+ k f; E - *)(X(rJ) - 

n=i KeT K 

772 

uEE E kk? up fl'(rr'JK •»**]. 

72—1 KeT ae£(K) 

<j=K\L 


(5.1) 


16 











for any 0 < r E C' 1 ([0, T] x 0), U E C' 1 ([0,T] x 0), U\gn = 0, where we have used notation \3. 1 71) for 
r n , U n and [KfMKW for U£, XJ% K , r v K , u£. 

We notice, comparing the terms in the “discrete” formula m with the terms in the “continuous” 
formula m, that Theorem o represents a discrete counterpart of the “continuous” relative energy 
inequality (E2D- The rest of this section is devoted to its proof. To this end, we shall follow the proof of 
the “continuous” relative energy inequality (see m and m and adapt it to the discrete case. 


Proof. First, noting that the numerical diffusion represented by terms (I4.2all~i~27i1i in the energy identity 
(14.11) is positive, we infer 

I\ + I 2 + I 3 < 0, (5-2) 

with 


a : = EJtHaki 2 - eFVrT). a ■■= 


KeT 


K&T 


I 3 := y^ (n f |V iT w n | 2 dx + (// + A) / |divu n | 2 dx). 
iftr Jk Jk ’ 


Next, we multiply the discrete continuity equation (I3.14al ) f- by b\^hK\ an< ^ sum over ^ € T to 
obtain 

h ■■= E YT {& ~ * _ ‘>|V»VI 2 = - E E ;kle? up [< • •vaK.kI 2 ■= A ( 5 . 3 ) 


KeT 


KgT o&£{K) 
a=K\L 


In the next step, taking — U n as test function in the discrete momentum equation (13.14bl) : we get 

h = - E { ^r(QK»K - e n K l ^K l ) • Uf. K = J 2 + J 3 + J 4 , 

KeT K 


with 


A=E E 

Ka'T ct€l£(K) 

<j=K\L 

J 3 = li ]T f Vu n : VUjf dx + (n + A) ]T f divu n dwU% dx 
A'eT K KeT K 

and 

A = - E E Ml 

KeT ae£(K) 
a=K\L 

We then multiolv the continuitv eauation (I3.14al)^ bv H' (rTT 1 ) and sum over all K E T and obtain 

- E - er'Itf'W- 1 ) = E E Me? up [< • n„, K ]H'(r" K -'). 

K&T K&Ta&£{K) 

a=K\L 


Observing that - q xi k 1 H'(rf < x ) = Qk{ h> {r\) - H' (r\ x )) + {o n K 

the last identity in the form 

h - - E nr(*«'«) - sr'^'W- 1 )) = j 5 + j 6 

KeT K 


cf K 1 )H'(rf < 1 ), we rewrite 


with J5 



E E Mer p K ■n„, K }H'(r' h r‘). 

K(fzT cr££(K) 
a=K\L 


(5.4) 
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Finally, thanks to the convexity of the function H , we have 

h ■= £ !§l[(rS - H[rl)) - (r^H'(r^) - H(r" h r 

K&T 

= £ - ipr”" 1 )) - £ - «■ - r^H’(r" K ~') - H( r”' 1 ) (5.5) 

AeT AeT 


k 

< E := *^7> 

AeT 

Now, we gather the expressions Q-& this is performed in several steps. 
Step 1 : Term _Zi + I4 + I 5 . We observe that 

n—1 1 j- 7-n—l |2 


I T T n I 2 n n \TT n \ 2 — n n ~ 1 \TT n ~ l \ 2 TT n ~ l + TT n 

]U h,K\ (jn _n— 1 \ ^aI^/i.aI ^A l U /i,A I , „ n -i U h,K + U h,K tTrn -l rrn a 

^ 1£>A - J - -2- H Sk -2-V°A,A “ U h,K b 

- - er^r 1 ) • u h< = ~(qk u k • k k - er^r 1 • - esr^sr 1 • ^ 


Consequently, 

t I T T — ^ IE 7„n |„.n rrn |2 _n—l|„.n—1 rrn—1|2 

-'l + 7 4 + -C - 2^ oWT" (£aI w A ~ “ Qk \ u K ~ U h,K I 

AeT^ " 

Tjn _ rr ™- 1 
1 u h,I< u h,K 


TT n 1 4- TP 1 
,C A,A + U /i,A _.n—1 


— U 


A 


- £ 1*1*- 
AeT 

Step 2 : Term Ji + J2. The contribution of the face cr = K\L to J\ reads 
jjn 4. f/ 77 _i_ 

-M Sk ’ 2 1 • (TV - W K • n„ lK ]+ - k| e ; '■• L ■ (C/Jl - TV [< • nV 

Similarly, the contribution of the face <r = K\L to J2 is 

■ (TV - UX,l)K ■ n,,K\* + M(?>! ■ (VJU - U£,k)K ■ *V,r] + - 


( 5 . 6 ) 


Consequently, 


J\ + h = - E E l a l£i 

Ke1~ (j=K\L(E:£k 


.nW ^MT + ^L 
V 2 


— u 


n, up 


Uh, aK • 


(5.7) 


Step 3: Term I 3 — J 3 . This term can be written in the form 

h ~ h = E [ l V -K' - U h )\ 2 dx + (m + A) [ |div(u" - U£) | 2 dx) 

AeT Jk Jk 

Ej (vc/h n : V (U% - w n ) + (/X + A) [ divttf div(ttf - u")). 
PTb 7 a v Jk > 


(5.8) 


AeT 


Step 4: Term ^2 + ^6 + T?. By virtue of (15.21) . (15.4115. 51) . we easily find that 


^2 + h + I7 — E (-^(^A I 7 ’a) — ^{Qk 1 I r a : )) ’ 


(5.9) 


AeT 


where the function i? is defined in (12.51) . 
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Step 5 : Term J5 + J§ + J7. Coming back to (l5.4H5.5p , we deduce that 


J5 + + J7 



«'«-')) + £ £ (5.10) 

K€lT ctG£(K) 
a=K\L 


Step 6: Conclusion 
According to (EZD-fESD, we have 


e —e 

2—1 2—1 


whence, writing this inequality by using expressions (15.611 (15.101) calculated in steps 1-5, we get 


E - *"‘K rl - v«T) + £ ' J Y{m<\r n K ) - Etervr 1 


i*i 


xer 


KeT 


+ E (m/ |V*(^-t/£)| 2 dz+( M + A) [ \dw(u n -UZ)\ 2 dx) 

KeT JK Jk 

< E (m f V x Uf :V x (Uf -u n )dx + (p + \) [ awu]}dw(uj;-u n )dx 

KeT V 7a 7a 


+ E M 

K&T 

- E E l^iei 

xeT <j=K\Le£ K 


TT n — TT n ~ l TT n ~ 1 4- T/ n 
-l U h,K u h,K / u hK + 


k 


n, up 




h,K ~ T '- J h,I< n -l 
— 


TT n -I- TT n 
u h,K + u h,L 


-u^]-U^ k [K-u^ k ] 


K I 


(5.11) 


- E E Ml«)[E4V • n,,*] + £ - Q n K )(H’{rl) - lT(r£- 

KgTc7=K\L&£ k KeT 

+ E E 

XeT tr=K\Lc£ K 

We obtain formula (ED by summing (15.111) ” from n = 1 to n = m and multiplying the resulting 
inequality by k. □ 


6 - Approximate relative energy inequality for the discrete problem 

The exact relative energy inequality as stated in Section 0 is a general inequality for the given numerical 
scheme, however it does not immediately provide a comparison of the approximate solution with the 
strong solution of the compressible Navier-Stokes equations. Its right hand side has to be conveniently 
transformed (modulo the possible appearance of residual terms vanishing as the space and time steps 
tend to 0) to provide such comparison tool via a Gronwall type argument. 

The goal of this section is to derive a version of the discrete relative energy inequality, still with 
arbitrary (sufficiently regular) test functions ( r,U ), that will be convenient for the comparison of the 
discrete solution with the strong solution. 

Lemma 6.1 (Approximate relative energy inequality). Suppose that P C R 3 is a bounded polyhedral 
domain and T its regular triangulation introduced in Section \3.11 Let the pressure p be a C 2 ( 0,oo) 
function satisfying hypotheses \1.5\) . 1 1. 61) with 7 > 3/2 and satisfying the additional condition 11.71) if 
7 < 2 . 

Let (g°,u°) € L^(fl) x W h (Q) and suppose that ( Q n )i< n <N € [L^(Q)] N , (u 11 )\< n < N € [W/^P)]^ is 
a solution of the discrete vroblem (13.141) with the viscositv coefficients a. A obevina E3P - 
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Then there exists 


c = c(M 0 , Eq, r, r, \p'\ C k \r,r], II (dtr, d(r, Vr, d t Vr, U, d t U, VE7, d t VU)\\ ioo^^sij) > 0 
(where r = max^ t x )^q^ r {t-, ®), £ = min, t x)eQf r (A x )^ such that for all m = 1,..., iV, we /iawe: 

£(g m , £/ m ) — <S(^°, u°|r(0), £/(0)) 

+ & E E (^ / \^x(u n - UJ ()\ 2 dx+ (fi + A) f |div(it n - Uf )\ 2 dx) 

n=l KgT JK JK 


<^EE W, V X UJ( : V x (Uf - u n ) dx+(p + A) / divUjfdiv(Ujf - u n ) dx 
n=lK&T 


IK 


IK 


_ n— 1 

— 1 h,K ^h,K 


TT n — TJ 


+‘EE l*l* 

n= 1 KgT 


+t££ £ H*“ p (t>E“ p -a; 

n= 1 A'gT ct££(K) 


TT n — i, n 
u h,I< ~ U K 


u; - us, K m:<- ■ n, x 


-*£ £ / P«)divt/“dx + *£ X / (r£-<&)^*W]“d* 


n=l A'gT 

m 


IK 


n=l A'gT' 


K 


' K 


a = ^ E £ l*l*' 

n=l A'gT 


Tjn _ TTn—l TT n ~ l 4- TT n 
-1 u h,K u h,I< ( u h,K^ u h,K n _i> 

o u k h 


n = -*EE E 


Fie 1 , 


n,up 


n=l A'eTo-=AT|Aef (AT) 


TT n + rr 1 
h, K + u h,L 


-&”™)-UZ iK [vZ-n a , K ] 


r 4 = - fc E E E 

71=1ATgT a=K\Le£(K) 

n = * £ E l*K* - *)—*’ “ ' ff ' Krl) 

n=l A'gT 


k 


( 6 . 1 ) 


/■ n" 

- A; E E / -fv\r n K )u n • Vr n dx +i^ fc +G m 

n=lKGT JK r K 

for any pair (r, [/) belonging to the class \3.18\) . where 

\G m \<ckJ2£(0\u n r n ,U n ), \R% k \<c(Vk + h A ), andA=\ (6-2) 

and where we have used notation \3.1 r A) for r n , U n and \3. 7HT&I) for UJJ, U^ ; r^-, u”. 

Proof. The right hand side of the relative energy inequality (EH) is a sum where 

m „ 

T 1 = E E (/* / . ^x U h ■■ V x (Uf - u n ) dx + (p + \) div Uf div (Uf - u n ) dx), 

n=i A'gT V Jk Jk 

(- 


n ~ k E E E Mer p *'(*TK''v/d- 

7i=l A'gT a=K\Le£(K) 

The term Ti will be kept as it is; all the other terms Tj will be transformed to a more convenient 
form, as described in the following steps. 
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Step 1 : Term T 2 . We have 


t 2 = t 2 ) i+* 2 j i, with t 2j1 = kJ2Y, \ R \e n K 

n=l KeT 


TT n — TT 


-l^h,K 


n—1 
h,K 


rn 

(uZk a " dR v =‘EE R u. 


n=l KeT 


where 


n ,K_m ^{ufa-u&y _\ K \ n - 1 ([u n -u n ~ i ] h , K ) 2 

n 2.1 2 -/' 


k ~ 2 0K k 

We may write by virtue of the first order Taylor formula applied to function 1 1—)• U(t, x), 


[U n - U 


n—l] 


\h,K 


k 


1 

\K\ 


IK 


T r j' tn 

-k L 


in — 1 


(9iU(z, x)dz 


- /i 


d.x 



1 f tn r / 



w\Sk 

k i [dtJJ{z) 

^ Jtn -1 

(x)dz 

h 

<\x 


where we have used the property (I9.2UI) of the projection onto the space 14(17). Therefore, thanks to 
the mass conservation (13.151) . we finally get 


|j#*i < - 1 "„ 2 

1 2,1 1 - 2 


I Tv | fc 114 U 11 ioo ( 0)T;£/ «x>(n ; R 3 )) ■ 


(6.3) 


Let us now decompose the term T 2) i as 


t 2 , 1 = r 2j2 + r 2>2 , withr 2j2 = 

n=l KgT 


TTn _ tt 
h IT u 1 


-1 ^h,K 


n—l 

h,K 


in 

([4V-<), and* 2>2 = fe£*2 )2 , 


n= 1 


where * 2 , 2 = 

KeT 


Tjn _ T TTl—l 

u h is ej 7 


1 


WJsT 




K* 1 - ua, K ) - z m sk 


TT n — TT n ~ l 


1 '- J h,K 


h,K 


KeT 


(6.4) 

■(ul-'-ul). 


By the same token as above, we may estimate the residual term as follows 


1*2, 2 I < k cM 0 \\d t u\\l x{0T . wl ^. R3) + cM 0 1/2 ( £ l^ler'K -1 - <| 2 ) 1/2 ||9tC/|lLoc(0, T ; i o O ( n ; ]R 3)), 

KeT 

where we have used the Holder inequality to treat the second term; whence, by virtue of estimate (14.2a[) . 

I* 2 , 2 1 < Vkc(M 0 ,E 0 , \\(d t U,d t \7U)\\ Lao(QT . R i 2) ). (6.5) 

Step 2 : Term T3. Employing the definition (13.131) of upwind quantities, we easily establish that 
*3 = *3,1 + *3,1) 

m m 

with T 3A =kJ2 T £ kler p (“”' Up -AT)' ,7 A“”'"wJr. % = lEE 8 W 

n=l KeT cr££(/C) n =1 ceSint 


and *3 ; ’f = |<r|^ 
Writing 


I TTn _ TTn |2 
\ u h,K u h,L I 


|[/n _ £jn |2 

[w" • n^ K } + + \a\ol— h ' L h,K [w" • n CTiL ] + , Vcr = A|L € £ int . 


Uh,K - u h,L = m Ik - W,* + UZ tff ~ UZ, L , <7 = K\L G T int , 


employing estimates (19.11) and (19.221) .-i and the continuity of the mean value 77” = of 77^* over 
faces <r, we infer by using the Taylor formula applied to function x U n (x), 


1*3,11 < h 2 c||V77||^ (QT . R9) |a|(^ + Ql)\K\, Vu = K\L G T int , 
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whence 


1^3,11 < h c||V£/| 


'A°°(Q t ;R 9 ) 

< h c(M 0 , Eq, ||VL/||£oo(q t; 


( £ £. h\o-\(gx + 8 l 

KeT a=K\Le£(K) 


5/6 r m 

k 


£(£ £ %i 

71—1 T CT(E:£(K) 


U 


n 16 


l/3-i 1/2 


( 6 . 6 ) 


provided 7 > 6/5, thanks to the discrete Holder inequality, the equivalence relation (13.211) . the equiva¬ 
lence of norms (19.351) and energy bounds listed in Corollary 14. 11 

Evidently, for each face a = K\L e £; nt , tr™ • + w™ • n CT) £ = 0; whence, finally 


,AT 


(6.7) 


m = E M«r p («r p -CT)-( u l k-u;)k■ 

n=l Kfz.'T cr(z£(K) 

Let us now decompose the term T \37 as 

m 

23,1 = ^3,2 + ^3,2, With 2?3^2 = k £ i?-3 2) 

n=l 
m 

n,2 = k y, E E M* up (tC p - ^r p ) ■ (tc - t4V)“” ,up - 

71=1 A'eT ae£(K) 

K 2 = E E m* w (C"; 7 -*r p ) - (t” - v*v)(«*; - cr p ) ■ «»,*■. 

K&T ae£(K) 

By virtue of discrete Holder’s inequality and the first order Taylor formula applied to function x 1 —)■ 
U n (x) in order to evaluate the difference U™ — Uft K , we get 

2 \ !/2 


and 


i$ )2 | < c||VU| 


L°°(Qt;R 9 ) 


71, Up 


-71,up Tj-n,up 

U CT U h,a 


(£ £ %iia 

KeT (j{z£(K) 


(£ £ h \m 

KeT (tg£{K) 

,/,270, (E E m 


n,up 170 


a\ 


KeT cr££(K) 


ul - C’ UP 


1\ l/<? 


where ^ + 7 /^ + 7 = 1, 70 = min{ 7 , 2} and 7 > 3/2. For the sum in the last term of the above product, 


we have 


£ £ %l 

KeT cr£.£(K) 


ui - u n j up * < c £ £ h \ 


<?\ i<-<r 


K£l~ a££(K) 


< c 


(E E 

A'eT ad£(K) 


A'eT 


) 2'yq —o 

< ch ^0 


V^/fijR 3 )’ 


where we have used the definition (13.131) . the Minkowski inequality and the interpolation inequalities 
()9. 181 fDTTTJl) . Now we can go back to the estimate of R 7 ^ 2 taking into account the upper bounds (14.91) . 
(|4.12H4.l3ft . in order to get 

17^3,21 < c(M 0 ,E 0 , IIVU|| L oo(g T;R 9)) (6.8) 

provided 7 > 3/2, where A is given in 
Finally, we rewrite term X 3 2 as 


T’S .2 — T 3 3 + i? 3 i3 , with i? 3 )3 — k £ III 


3j 


n= 1 


UL 

7 - 3,3 = 1 E E E m*“ p K 7 -*r p ) • (a; - <av)£4”;7 • 

71=1 A'eT (tg£(K) 


(6.9) 


2?3,3 = £ £ Me; 

A'eT o-e£(A') 


)/ |- A r71,U 


,ti, up / rr n ’ u P _ ,^71,up 


T"-CT-)( 


-71,Up 1-7-71,Up 
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whence 


( 6 . 10 ) 


1 ^3,s| < c(|| VL/|| Loo(QtiR9) ) k £ S(g n , u n | r n , C7"). 


n=l 


Step 3: Term T 4 . Using the Stokes formula and the property (19.231) in Lemma 19.21 we easily see that 


T 4 = -kJ2 / p(6k) div un dx 


n=l K&T' 


K 


( 6 . 11 ) 


Step 4: Term T 5 . Using the Taylor formula, we get 


H’m - H'ir^ 1 ) = H"(r n K )(r n K - r^ 1 ) - ~-H'"{r n K ){r n K - r^ 1 ) 


■w /—n \ / n n—l\2 


where r 1 ]^ € [min (r v K , r^),max(r^- , r^)]; we infer 


T 5 = T 5)1 + with T 5>1 = kY,J2 \ R \( r K - e n K) F "n ,,K k K > ^5,1 = ^E ^ 5 ,f, and 


rJ(r n \ r n r " _1 

p Ua'1 r K - r I< 


Rif = f|A-|g'"(n.) VK ' (<& - rfc). 

Consequently, by the first order Taylor formula applied to function t eT r(t, x) on the interval (t n _i, 
and thanks to the mass conservation (13.151) 


n=i A'eT 
„n „n—1\2 


'a: 


n=l A'eT 


( r n _ 

—n \ v AT r A 1 


l-R5.1l < k c(M 0 ,r,r, Ip'bqfcr], l|du , ||z,°°(Q T )), 

where r, f are defined in (13.181) . 

Let us now decompose T^i as follows: 

p'ir k) 1 


( 6 . 12 ) 


n, 1 = r 5 , 2 + T? 5 , 2 , with T 5 , 2 =1’EE / (*$ - ^^r^ufdx, i ? 5 , 2 = £ E Rs,f, and 


-^5,2 ~ 


n=l A'eT 


/a: 


’A' 


n=l A'eT 


(6.13) 

In accordance with (13.171) . here and in the sequel, [ dtr] n (x ) = dtr(t n ,x). We write using twice the 
Taylor formula in the integral form and the Fubini theorem, 


I r>n,K< _ 1 
l R 5,2 I - 1 




P\r n K )r n M - rl) I I / dfr(z)dzdsdx 
J J t n — 1 J S 


< 
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^ In 11^ - r I \LHK) 

'K 1 

Therefore, by virtue of Corollary 14.11 we have estimate 


9fM-)llA7'(X) dzds - 


IR 5 . 2 I < k c(M 0 ,E 0 ,r,r, |p'| c i([rT]’ ll 9 t r ILi(o,T ; LV(n)) 


(6.14) 
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Step 5: Term Tq. Using the same argumentation as in formula (16.71) . we may write 


m 

T 6 = T 6)1 + Rs tl , R 6 ,i = kY,Y, E R gJ’ K ’ with 

n=lKeT aeS{K) 

m 

T 6,i = k E E E M 6k - ff'(r£ _1 ))«S • n a>K , and 

n =1 K&T a=K\Le£(K) 

jyn./T. K 
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ryn,a,K _ |„.|/_ra,up 
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=K\Le£(K) 
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(6.15) 


We estimate this term separately for 7 < 2 and 7 > 2. If 7 < 2, motivated by Lemma 14.21 we may write 

I T)n,a,K I 
^1 


f^r Jir |<V^||Vfr , (r)|| Loo(gT;R 3 ) |a| 

|~n,up _ n n I ,- ,- 
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^a,K | 


+ |^’ UP - 2aI ' n ^l)> ( 6 ' 16 ) 

where we again use the Hrst order Taylor formula applied to function H' 
and where the numbers are dehned in Lemma 14.11 on 

Young inequalities yields 


ied to function H' between endpoints r\ 1 , r™ 1 , 
Consequently, an application of the Holder and 

m _ _ („n,up _ _n \2 \/2 
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Ve deduce employing the discrete Holder inequality 
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}E (E l*l sf-^f' 3 ] 1/2 [i- E (E MM-T) V3 ] V2 - 
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' back to (TB7 


1^6,i| < \/hc(M 0 ,U 0 ,r,r, |p'| c . ([riF]) ,||Vr|| Loo(QT;R 3 ) ) 


(6.18) 
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provided 7 > 12/11, where we use estimate (14.161) . estimates (14.101) . (14.121) of Corollary 14.11 and equiva¬ 
lence relation (19.351) . In the case 7 > 2, the same final bound may be obtained by a similar argument, 
replacing the estimate (14.161) by (14.151) . 

Let us now decompose the term Tg ) 1 as 


n, 1 = r 6 , 2 + r 6 , 2 , with t 6 , 2 = k £ £ £ - C _1 )K • 

n=lKeT <r=K\Le£(K) 

m 

Re,2 = & E E E and 

n=l KelCa££{K) 


R’St K = MeS/ff'MT 1 ) - ff'W 1 ) - ff"W- 1 )W 1 - ’T 1 )) 


U <7 • 


Therefore, by virtue of the second order Taylor formula applied to function H’, Holder’s inequality, 
(|9.30l) . (19.351) . and (14.91) . (|4. 121 ) in Corollary 14.11 we have, provided 7 > 6/5, 


|-R6,2| < hc(\H \c([r,r]) + \H |c([r,r])) ll^ r llL 00 (Q T ;R 3 ) lll ) lli° o (0,T;L'i'(n)) ll' u llL 2 (0,r;V2(Q;R 3 )) 

<hc(M 0 ,E 0 ,r,r,\p'\ c i(\r,r}), l|Vr|| La c ( Q T;H 3j), (6.19) 

where in the first line we have used notation (13.111) . 

Let us now deal with the term T@ )2 . Noting that / Vr ™^ 1 dx = E \ a \( r a 1 ~ r K 1 ) n a, K , we may 

K <t££(K) 

write 


E - T-r 1 )K • n,,*] 

a££(K) 


= -[ g n K H"(r n K - 1 )u n K -Vr n - 1 dx+ £ “ C“ X )« - O ■ 

JK 

Consequently, Tg )2 = + i?6,3 5 with 


ae£(K) 


r 6,3 = ~k E E L • Vr - 1 dx, 

n=l K&T" " 


IK 


Re,3 = k E E I,. Q n KH"(r n K - l ){u n - <) • Vr^ 1 dx 
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IK 


71=1 K(zlT ct(E:£(K) 


|i? 6 , 3 | < c||iL"(r)Vr| 
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*EEI 

71=1 KeT 


\Qk\ 
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fc E E E HeSrlU'«>(«’)ll w S- w Srll L 7j (A -. R 3 ) J> 

n=lK£Ta££{K) 


70 = min{ 7 , 2 }, 


where we have used the Holder inequality, and also the Taylor formula applied to function x i->- r(t n _i, x) 
together with equivalence relation (13.211) yielding \a\h < \K\, to treat the second term. Consequently, 
by virtue of Holder’s inequality, interpolation inequality (19.51) (to estimate ||u n — u^-H^/ k3 by 
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/j( 57 o 6)/(270)||V x u n || L 2( A -. R 9), 70 = min{7,2}) in the first term, and by the the Holder inequality and 
(I9.5H9.6I) (to estimate ||u" — RkIIlT (K- 9 3 ) by /j( 5 'to- 6 )/( 2 7 o)Hv^u” \\l 2 (k-m. 9 )) i n the second term, we get 

TLl /2 ^ /2 

< ch^-o >/(»») l|ff"(r)Vr|| 1 „ (gi . iR , ) fc £ ( E ll*ll!™ (K) ) '( E IIV.ulV;*.)) 

n =1 KgT KgT 


< c/i (5 ' T °- 6 )/ (2 ' 7 °) ||if"(r)Vr| 


RE 


( X! WSKW'L'yo(K) 


l/7o 


(E ilv 


a II L 2 (K;R 9 ) 


1/2 


n= 1 Ker -ffeT 

provided 7 > 6/5, where we have used the discrete Holder inequality and the algebraic inequality (19.71) . 


Now it remains to use (14.91) . (14.121) in Corollary 14. II in order to get 

1^6,31 < h A c(M 0 ,E 0 ,r,r, |p , |c 1 (fc.rl)ll Vr lli < » ( gT ; R3 ) ), 

where A is defined in (16.21) . 

Finally we write = Tq ,4 + Re, 4 , with 

1 m = -*£ £ f dz, 

v^1K€T Jk r K 


n- 

m 


(6.20) 


#6,4 = k J2 £ /. e n K(H"{r^)Vr n - H"r^ 1 ) ■ u n dx, 
n=lK£T JK 

where by the same token as in (16.141) . 

\R 6A \ < k c(M 0 ,E 0 ,r,r, \p'\ c ^[r,r]), ||'Vr, II^Vr|| i2(0iT;i67/( 5 7 -6) (n . R 3 )) ). 


( 6 . 21 ) 


( 6 . 22 ) 


We are now in position to conclude the proof of Lemma 16.11 we obtain the inequality dSUD by 
gathering the principal terms (16.41) , (16.91) , (16.111) , (16.161) , (16.211) and the residual terms estimated in (16.61) , 
(16.51) . (16.61) . (16.81) . (16.101) . (16.121) . (16.141) . (16.161) . (16.171) . (I6.19|) , (16.201) . (16.221) at the right hand side 

J2i=\Ti of the discrete relative energy inequality (15.11) . 

□ 


7 - A discrete identity satisfied by the strong solution 


This section is devoted to the proof of a discrete identity satisfied by any strong solution. This identity 
is stated in Lemma O below. It will be used in combination with the approximate relative energy 
inequality stated in Lemma 16.11 to deduce the convenient form of the relative energy inequality verified 
by any function being a strong solution to the compressible Navier-Stokes system. This last step is 
performed in the next section. 

Lemma 7.1 (A discrete identity for strong solutions). Suppose that O C l 3 is a bounded polyhedral 
domain and T a regular triangulation introduced in Section HOI Let the pressure p be a C 2 ( 0,oo) 
function satisfying hypotheses and m with 7 > 3/2. Let (r, U) belong to the class \3.18\) satisfy 
equation «nn with the viscosity coefficients fi, A obeying P^\ )- 

Let (g°,u°) € L^(fi) x W h (Q) and suppose that ( Q n )i< n <N € [L^(Q)] N , (u n )\< n <N € [W h (Q)] N is 
a solution of the discrete problem (16.141) . Then there exists 


c = 


c(m 0 , E 0 ,r,r, bVfcr]), II(Vr, d t r, U, VE7, V 2 U, d t U, d 2 U, d t VU, )|| L oc (Qt;R 58 ) )) > 0, 
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such that for any m = 1,..., N, the following identity holds: 


m ,, 

k E E / u h • V(«" - m + (/X + A) div L7? div(u" - 17?)) dx 

n=lK£T JK 


+ k H E 


/A" 


Tjn _ TJ n ~ l 

rp KK k KK dx 


n=l KeT" 
m 

+ k Y.Y. E MC" P K'7■ n«,K](tc-t^,jr)• («’“ p - 

n=lKGT(7G5(^) 

m m „ 

+ ^EE / i«) div IT” dx + fc£ E / pVi^'V/Mx+ft^O, 


(7.1) 


n=l ATeT 


'A' 


n=l A'eT 


'A' 


where 

< c(/i + fc) 

and where we have used notation 4,7.7 71) /or r n , U n and 4,7.71I V. <7I) for U)j, U’/ A -, rf^, u”. 

Before starting the proof we recall an auxiliary algebraic inequality whose straightforward proof is 
left to the reader, and introduce some notations. 

Lemma 7.2. Lef p satisfies assumptions (ESI) and (ESP- Lei 0 < a < b < oo. Then there exists 
c = c(a , b) > 0 such that for all g £ [0, oo) and r € [a, 6] f/iere holds 


E(g\r) > c(a,&)(l fl+ \[ a / 2 , 26 ]( 0 ) + £ 7 l,R + \[a/ 2 , 2 b] ( 0 ) + ( 0 - r) 2 l[a/ 2 , 2 b](£) 


where E(g\r ) is defined in \2.5\) . 

If we take in Lemma 17.21 £ = g n (x), g n € L^(fl), r = f n (x), a = r, b = r (where r is a function 
belonging to class (13.181) and r, r are its lower and upper bounds, respectively), we obtain 


L;(^(x)|r-(.T))> C (r,r)(l R+ \ fc/2 , 2F] (^(x))+(^)^(x)l R+V[ , /2 , 2?] (^(.T))+(^(x)-r"(x)) 2 l fc/2i2F] (^(x))) 

( 7- 2) 

Now, for fixed numbers r and r and fixed functions g n € Lt (fl), n = 0,..., TV, we introduce the residual 
and essential subsets of fl (relative to g n ) as follows: 


Kss = {xcn E < g n (x) < 2If}, Nfi es = to \N™ 


(7.3) 


and we set 

bless (x) = g(x)l NSss (x), hires (x) = g(x)l N n s (x), x € fi, g^L l (yi). 
Integrating inequality (17.21) we deduce 




dx < £(g n , u n 


r n , U n ). 


(7.4) 


for any pair (r, U ) belonging to the class (13.181) and any € Lh{ 12). 

We are now ready to proceed to the proof of Lemma 17.11 

Proof. We start by projecting the momentum equation to the discrete spaces. Since (r, U) satisfies (11.11) 
and belongs to the class (13.181) . Equation (ll.lbl) can be rewritten in the form 


rdtU + rU ■ V£7 + Vp(r) = gAU + (g + A)V div U. 


(7.5) 
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We write equation (17.51) at t = t n , multiply scalarly by u n — Ufi, and integrate over fL We get, after 
summation from n = 1 to m, 


b m . 

^71 = 0, with Ti = -kY J / (//A[/ n +(> +A)VdivC/ n ) ■ (u n - £/£) dx, 

i= 1 n= 1 n 

m „ m ,, 

r 2 = ky / r n [d t u] n • (' u n - U£) dx, r 3 = ky / r n [/ n • Vt/ n • (w n - U%) dx 

^l Jn n^l Jn 

m « m 

T 4 = ky / Vp(r n ) • u n dx, T 5 = -ky / Vp(r n ) • U% dx. 


n=l 1 


71=1 ’ 


(7.6) 


In the steps below, we deal with each of the terms 7). 

Step 1 : Term 7). Integrating by parts, we get: 

71 = 71,i + 7^1,1, 

771 r. 

with 71,i = kJ2 y / (/zVl/£ : V(w n - U%) + (ji + A) div E7£ div(it n - t/£)) dx 

n=LK£T JK 
and 72-1,1 = ii + 72, with 

m ,, 

h = ky E / {/xV(?7--C/^):V( W "-I7^) + ( M + A)div((C/"-C/^)div(n"-?7^))dx ( 7 . 7 ) 

n =lKeT JK 


h = -fc E E E / • VC/ n • («" - U h) + (A + M)divt/ n (u n - UZ) ■ n a y d5 

n=lKeTcre£(K) Ja 


rn « 

-*E E / • VET 


71=1 crE£ 1 


- C/T 


+ (A + //) div U r ‘ 


w n - t/T 


■ n<r) dS, 


where in the last line n CT is a unit normal to er and n<r is the jump over sigma (with respect to n CT ) 
defined in Lemma 19.51 Employing estimate (19.221) we easily verify that 

|/i| < h c(Mo,Eq, ||VU|U=o (0iT;W ,i,oo (n)) ). 

Since the integral over any face u of the jump of a function from V/i(fi) is zero, we may write 


h = ~kyy \JJ-n 0 

n= 1 ae£ Ja 


(yu n - [vu r 


u n - UH 


J cr, no- 


+(X + fi)( div U n — [div U n ] ( 


u n - un 


n a ) clS 1 ; 


whence by using the first order Taylor formula applied to functions x >-)• VU n (x) to evaluate the 
differences \7U n — [VU n ] a , divU 71 — [div[/ n ] ff , and Holder’s inequality, 


\I 2 \ < kh c ||V 2 LI|| L oo ( q t;K 27 ) y y \J\a\\Zh (~^= 


71=1 (7G£ 


1 


u n -un 


L 2 (ct;R 3 )> 


771 J 

< kh c||V 2 ?7|| i ^ (QT . R 27 ) ^ (M h + - 


71=1 <7G£ 


u n -un 


2 

L 2 (o 


Therefore, 


|72i,i| < hc(M 0 , E 0 , ||U, VE7, V L/|| LO o(q T) 
where we have employed Lemma 19.51 and (14.91) in Corollary 14.11 


(7.8) 
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Step 2 : Term 75- Let us now decompose the term 75 as 


75 = 72,1 + 'R-2,u 

m 

with 75,i = & E E 

n=\K&T - 


m—1 


m 


K 


Tjn _ Tjn 

r n ~ l - t - (u n - UJi) dx, 77 2 ,i =^E Kf, 


and 77”f = J^(r n - r n ~ l )[d t \l] n • ( u n - U%) dx + r n ~ l {[d t U] 
The remainder can be rewritten as follows 


n= 1 KgT 
jjn _ jjn-l 

” - - k - )•(«*”- {/?) dx. 


/-j-j tl^K 

= 


K Jtn -1 


r(t,-)dt [d t \]] n ■ (u n - U£) dx 


+ k] K r 


,71 — 1 


‘'n / 


Jtn — 1 J & 


dfU(z, -)d2;ds • («" — U%) dx; 


whence, 


17^2’fl < k (|HIl°°(Q t ) + ||^r|| i oo ( Q T ))(||9tL/'|| £; oo(Q T . R 3)|iL| J / 6 (||w' t || Z/ 6( A ') + \\U%\\ L 6( K )) 

+ ||0|[/ ra || L 6/5(n ;K 3))(||u n || L 6(X) + re|| L 6 (J0 ). 

Consequently, by the same token as in (16.141) or (16.221) . 

|77 2 ,i| <fcc(Mo,^,f,||(a t r,t7,a t t7,Vt7,)|| Loo(gT;R i6 ) ,||^t7 llL 2 (0,T;L 6 / 5 (n ; R 3 )))> (7-9) 

where we have used the Holder and Young inequalities, the estimates (111.211) . (19.241) . (19.291) . (19.MOD , and 
to the energy bound (14.91) from Corollary 14.11 

Step 2a: Term 75,i- We decompose the term 75,1 as 
75,1 = 75,2 + 7^2,2, 

171 r TT n _ T J n ~ l 171 

n—1 ^ ^ /_ _7i r tti\ j __ ^r» 7 \ A \ ' ^r, n,K 


with 75,2 = k 

n=l KgT' 


K 


' K r 


and 77o’<f : = / (r n 1 — r 


71 — 1 „7l — 1 


£/ n - U 


n —1 


-2,2 


IK 


K 


• K - UJt) dx, 77 2 ,2 = kJ2J2 n V2 , 

n=lKeT 

(u n -U%) dx; 


therefore, 


177-2,21 — I E ^2,2 I — ^ C ll^ r llL°°(Q T ;R 3 )ll^f7|lL°°(Q T ;]f 

K&T 


\n n — TT n II 

l W u h IIl 6 (0;IE 


Consequently, by virtue of formula (I4.10|) in Corollary 14.11 and estimates (I9.3U|) , (19.251) , 

|77 2 , 2 | < hc(M 0 ,E 0l ||(Vr, 7/,9 i [/,VJ7)|| L oc (QT . R i8 ) ). 

Step 2b: Term 75,2- We decompose the term 75,2 as 
75,2 = 75,3 + 772,3, 


(7.10) 


with 75,3 = & E E 

n =1 ifeT 


r 


77—1 

^h,K 


Tjn _ Tf 
^ h IT '-J i 


K 


and TZjJ = 


+ 




K " A; 

V n_ jjn-l U n_ U n-l. 


(u n - Uh) dx, 77 2 ,3 = ^ E E ^2,3 

n=lK&T 


k 


k 


- h 


u n - U'O) dx 



s 

1 

C3 

3 

1 


lU n - u n ~ 1 ^ 

) 

k 

h 

k 

h,K 7 


) • (w n - U%) dx= /f + /f. 
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We have 


I t k I _ 1 

1^2 I - T J K 


r /r rt n -i r r^n >. 

/ ( / dtU(z)dz — / <9t£/(~)d-3 )-(u n — Uh)dx 

Jk VL Ji n _i L A n _i 


< E 71-1 

~ k rR . 


in— 1 


V, 




- h 


L 6 / 5 (K;R 3 ) 


Ufr H£6(Jf; R 3 ) 


where we have used the Fubini theorem, Holder’s inequality and (19.11) . (|9.22|) „—i . Further, employing 
the Sobolev inequality on Crouzeix-Raviart space (19.HOI) . Holder’s and Young’s inequalities and estimate 
(19.221) 8= i, we get 


Ji r 

E \%\ < tJr- 1 ll«" - ^IIl6 (0;R 3 ) / V^L/(z) 

k'cT h Jin-l 


KeT 

< 


L 6 / 5 (C;R 3 ) 


dz 




+ 


h r tn 2 


• 3 )^ k 


in— 1 


y x d t U(z) 


L 6 / 5 (n ; R 3 ) 


We reserve the similar treatment to the term . Resuming these calculations we get by using Corollary 

ro 


|7*2,31 < h c(M 0 , A(), || (n £7, VR, t/) II L°°(Q t; R 16 ), ||StV£/|| L 2 ( 0 ,T;L 6 / 5 (n ; R 9 )))- 
Step 2c: Term 72 3 - We rewrite this term in the form 


(7.11) 


72,3 = 72,4 + 77-2,4, 77-2,4 = ^ E 77-2,4 , 

n=l KeT 


with 72,4 = k EE 


_ n—1 
—1 '~ J h,K ^h,K 


TT n — TT 


n= 1 KeT 


IK 


1 K 


TTL — 1 


/. rr« _ TT r 

and 1Z n ’ K ~ f r n ~ 1 ° hJ< Uh ’ K 

and 77-2 4 ~ J k t k k 


( U K - Uh,K ) dx, 


(<*•” - «*) - (w - Kk)) 


First we write, as in 


[U n - U n ~ v 


h.K 


E / f- 

AT 7k Lfc L 


'in— 1 


dtU(z, x)d z 


- h 


dx 


|A"| Jk '-k 


’tn — 1 


[5tU(^)j ^(x)dzjdx < ||[<9 t U] /i J| L oc (0 , T;i oo (Q;R 3 )) < ||5tU 


(7.12) 


L oo (0,T;L oo (S2;R 3 )) 


Next we evaluate u n — u^ employing (IQ.llh -o. and Uf t — Ufi K by using (19.11b ,— (19.221) s =i. Finally we 
employ the Holder inequality to get 

17*2,41 < h c(Mq, Eq, r, ||(c?tL7, U, VE/)|| Loc ,(q t . r i 5), ||RtVL7)|| i 2(o, t ; l 6 / 5 (0;R 9 )))- (7-13) 


Step 3: Term T?,- Let us hrst decompose 7jj as 
75 = 73,i + 77-3,1, 


with r 3 ,i = k E E / r K u h,K • vc/ n •« - uIk) dx, 7e 3 ,i = k E E Kf 


77=1 KeT 


IK 


77=1 KeT 


and TT.o’f' = f {r n - r£)[/ n • VC/ n • (u n - U?) dx + / r£(t/ n - L7?) • VC/ n • (u n - U?) dx 
’ 7k 7k 

+ j K r n K(UZ - U^ K ) ■ VE7" ■ (u n - UJi) dx + f K r\Ul K ■ X7U n ■ (u n -u\- (C7 1 - U% K )) dx. 
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We find that 


\K. 


n,K | 
3,1 I 


< h 


\K\ 1/2 ( II 


u 


IIl 2 (K;R 3 ) + IIl 2 (K;R 3 )) + |-h. | 1/2 (|| Vu" IIl 2 (AT;R 3 ) + \\^ U h\\ L 2 (K;S. 3 ))] 

2 

X (IMIl^Qt) + ll^ 7 r llL°°(Q T ;R 3 )) (|I^IIl°°(Q t ;R 3 ) + II^ 7 ^ 7 "IIl°°(Q t ;R 9 )) > 


where we have used several times Holder’s inequality and the standard first order Taylor formula ( to 
evaluate r n — r^), along with the estimates (19.211) (to evaluate U n — Ufi), (19.11) . (19.221) „_i (to evaluate 
U% — U% K ), (19.11) (to evaluate u n — u v K ). 

Consequently, using again (19.221) „-i (to estimate ||V£/^|| L 2(^. R 3)), the definition of | • \v£(£i) norm, 
the Sobolev inequality (19.1101) and the energy bound (14.91) from Corollary 14.11 we conclude that 


1^3,i| < h c(M 0 ,E 0 ,r, ||(Vr, U, VU)\\ L , 


(Qt;K 15 ))' 


Now we shall deal wit term 731 . Integrating by parts, we get: 


(7.14) 



E wicPiK ■ n«,K\u: ■ w-ufa) 

ae£(K) 


E W\r n K [Ufa-n a , K ](UfaUfa)-(u n K 

ae£{K) 


Ufa) 


thanks to the the fact that J2ae£(K) la Ufa ' = 0. 

Next we write 


73,1 — 73,2 + 7^3,2, 7^-3,2 — k E 7^3,2; 

71— 1 


r 3 ,2 = k J2 £ £ MC’ up K', up ■ n„, K ](t/; - i4» • (*;« - £C P ), ( 7 . 15 ) 

7i=l KeT a^£(K) 


and 7^3 2 = E E MMr - ' n.,A'](t/ ff n - Ufa) • « - Ufa) 

Kfz.'T ae£(K) 

+ £ £ MC up [(HU - £>E up ) • n„, K ] (t/; - u^ K ) ■ « - t/; iK ) 

KaT cr(z£(K ) 

+ E E MC“ P (CT' n^K<A” - *AV) ' (Kr - <;7) - (VI.K - <D) • 

Ke 7~ a^£(K) 

We may use several times the Taylor formula (in order to estimate r^ — f” ,up , U™ — UJJ K , UjJ K — U^’^ p ) 
to get the bound 

1 7 ^ 3 , 2 ! < hc\\r\\ w i,oo^(l + l|t7’llw 1 .°°(n;R 3 )) E ^MI u a'I 


KeT 


n - u n I 

U (T 5 


/ \ Z 

+c||r|| 

+ ll^llw/O.^dniR 3 ) ) E E %H< 

KeT aE£(K) 

where by virtue of Holder’s inequality, (19.101) . (19.311) . (19.181) (19.191) . 

E h M\ u x\ <c( E %IKI 6 ) 1/e < c[( E ll“" - “*:ll V*»)) I/e 

KeT aeT K&T 
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+ 


+ 

(E 

u n 1 

6 5 1 / 6 ' 

L 6 (A; R 3 ) ) 


AeT 



E 

E 

h\a 

II-ji 71 i, n \ < 

\\ U K ~ u a\ E 

AeT a e£(A) 



E 

E II 

u n - 

7/ n || 2 

“o-Hl 2 (A;R 3 ) 


< c 


(El 

KeT 


I V U r , 


IL 2 (A' ;I 


1/2 


1/2 


KeT 

, V 2 i , 

< he 


KSL'T cr££(K) 


(E 

KeT 


I Vu 


n llL 2 (A;R 9 ) 


1/2 


Consequently, we may use (14.91) to conclude 

1 ^ 3 , 2 | < hc(^M 0 ,E 0 ,r, ||Vr, C7, V17 || LOO (q t;R i5 ) ). 


(7.16) 


Step 4: Terms Ta and T 5 . We decompose Ta as 


Ta — Taa + 77-4, 1 , 

m r. m » 

with Tip = k E E j p'(r r ^)u n ■ Vr n dx, 77.44 = k E E J (p'(r n ) — p'(r r ^)^ju n ■ Vr n dx; 

(7.17) 


n=l AeT 


n=l AeT 


IK 


whence 


17*4,11 < hc(M 0 ,E 0 ,r,r, |pV([r,r]), ll Vr llL°°(Q T ;R 3 ))- 
Employing integration by parts, we infer 

m r. 

% = r 5) i + 77-5^, with r 5) i = -fc E E / v ^( r?l ) ■ d * 

n=lKT JK 

m « 

7*5,1 = ^E E_4v P (o• (^ n - uji)dx 


n =1 AeT 


and 


172 - 5 , 1 1 < ^c(r,|p , | CM ,||Vr,Vl/|| Loo(QT;R i 2 ) ). 
Integrating by parts, we obtain 

m ,, 

r 5 ,i = fc E E l,p{r n )divU n dx- 

ri=\ KeT’ 


IK 


whence 


75,1 = 75,2 + 77-5,2, with 75,2 = fc E E l T _p(r?<)divU n dx, 

n=l A’eT" ~ 


/ft" 


77.5,2 = -fc E E LM r K ) - p( r?l )) div d3; 

n=l KdiT" ' 


IK 


and 


(7.18) 


(7.19) 


(7.20) 


172 - 5 , 2 1 </ic(r,|p , | fci? ],||Vr,Vl7|| LO o ( Q T;R ia ) ). (7.21) 

Gathering the formulae (17.71) , (17.121) , (17.151) , (17.171) , (17.191) and estimates for the residual terms (17.81) , 
(17.91 17.13D , (17.14fT7TT5T) , (17.181) , (|7.20l) , (17.211) concludes the proof of Lemma 17.11 □ 
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8 - End of the proof of the error estimate (Theorem 13.11 ) 


In this Section we put together the relative energy inequality (ED and the identity ED derived in the 
previous section. The final inequality resulting from this manipulation is formulated in the following 
lemma. 


Lemma 8.1. Under assumptions of Theorem \3.1\ there exists a positive number 

c = c(M 0 ,E 0 ,r,r, \p'\cp\r,r]), ||(Vr, d t r, d t Vr, d?r, U, VU, V 2 t/, d t U, $VE/)|| L oo ( q t . R 85 ) 
(depending tacitly also on T, 9 q, 7 , diam(O), |0|j, such that for all m = 1,..., N, there holds: 


£(g m ,u m \r m ,u m )+k^J2J2 / IV x (u n - U£)| 2 dx 


n =1 KeT 


IK 


< 


h A + Vk + £(g°, 


u° |r°, 


U°) +ckY J £{Q 

n =1 


W 


where A is defined in IQII . 

Proof. Gathering the formulae ED an( i ED i one gets 


£{g m , u m r m , U m ) - £(g°, u° r(0), 17(0)) + gk £ W ~ U h 


n =1 


V 2 (Q;R 3 ) 


< V\ + V 2 + 7*3 + Q 


( 8 . 1 ) 


where 


= E i^Ker'-'x 

n =1 KgT 


TT n — 77 n_1 
-U U h,K u h,I< 


(u, 


h,K — U> K I ) 


*>-*E E E ki(er p -c up )(GT p -< ,u dd a "- u ‘V)GT' n -'. 

n=l ifeT a=K\Le£(K) 


K, 


m „ 

v 3 = k E E /. (p(dr) - Kf£)) div ^ d * 

n=l KgT^ K 


1=1 KeT 

m 

+ k T. E [j 

n=lK&T 

Q = n% k + R% k + G m . 


r n „n 

1 1< 9 k J 


K r n K 


p (r r fi)u n ■ V? ,n dx + 


' K 9k I 


'k rl 


p \r‘j C )[dtr] n dx 


Now, we estimate conveniently the terms V\, V 2 , V 3 in three steps. 
Step 1 : Term V \. We have 


TT n — TT 
u h,K u ‘ 


n —1 
h,K 


< 


' tn — 1 




L°°(A';R 9 ) dz — 


||3t?7(z)|| L oo( A ';R 9 ) d Z, 


where we have used (19.201) . 
According to Lemma 17.21 

with any p > 1. In particular, 


6 - A ll R+\\r/2,2r\{Q) < c(p)£ P (£>|0 

le> - ^I 6/5 i« + \|v/ 2,25=](^) < c£(d r ) 


( 8 . 2 ) 
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provided 7 > 6/5. 

We get by using the Holder inequality, 


E \ k \(Qk 1 ' K 

KeT 


_ n —1 

I l 


Tjn _ TT 

u h u u ■> 


(Uh,K ~ U K ) E c (|I^^IIl°°(Q t ;R 3 ) + II^^ 7 ^IIl°°(Q t ;R 9 )) X 


(£ iKiier 1 -’-r 1 i 2 ifcA2fife)) 1/2 + ( e wi!*' 1 - ’•r 1 i 6/5 iR + \fc/2,2fi(w)) 


5/6 


ATeT 


(E 

A'eT 


A'eT 


\K\ 


TT n — n n 
u h,K U K 


6\ 1/6 


< 


c(||(a t J7,a t V(7)|| i oo ( Q T;R i2 ) )(£: 1/2 (^ n ~ 1 ,u n - 1 |r n - 1 ,J7 n - 1 ) 


+f 5 / 6 (£> n_1 ,ti n_1 |r n_1 , £/ n_1 )) ( E ||C/£ 


A'eT 


— ?/ n II 6 

A' W aIIl6(A';R 3 ) 


\ 1/6 

J ’ 


where we have used (|8.2D and estimate (|4.13l) to obtain the last line. 
Now, by the Minkowski inequality, 


( 51 W U h,K ~ u k\\%( K -,R 3 ) 


k&t 


1/6 


< 


(E \m, K - ui) - (u% - u n )\\% {K]i 


+ 11 U?-u n ' 


K&T 


-TO 3 'i E c 


1/6 


lA6(0;R 3 ) ^ 


U n - u? 


y h 2 (n ; R 3 ) 


where we have used estimate (19.141) and the Sobolev inequality (19.301) . Finally, employing Young’s 
inequality, and estimate (14.131) . we arrive at 


\Vi\ < c(S,M 0 ,E 0 ,r,r,\\(U,VU,d t U,d t VU)\\ Loo{QT ^ 5) ) 

m m 

x (kS(e°, U°\r°, U°) + k E £(e n ,u n \r n , U n )) + 5k E \™ n ~ U r h 


n= 1 


n=1 


v. 2 (n ; R 3 ) 


(8.3) 


with any <5 > 0 . 

Step 2: Term V 2 • We write V 2 = k E™=i ^2 where Lemma 17.21 and the Holder inequality yield, 
similarly as in the previous step, 


1^2 I E c(r,r, ||VL/"|| i0 o(Q T . R 9))x 

E E m h(E i ^' up ,c up ) + ^ 2/3 (^’ up ,c up ) |t>E up i i^T - ^ ,up i 

A'eT ctg£(K) 


< 


+ 


c(r,r,||(t7,VC/)|| Loo{ Q T . R i 2) )[( E E W\h(E(g™\C' ip )) 

KeT a€.£(K) 

( Y, E MhEieyiC^)) 213 

KaT cr^£(K) 


(E E \ a \ h 

A €7” a(z£(K) 


1/2 


T jn,up ~ n , up 

U (J 


6s 1/6 


provided 7 > 3/2. Next, we observe that the contribution of the face a = I\\L to the sums Ea reT 
E CTe£W |a|/ l £?(er p l^’ up ) and Exer T,ae£(K) W\h\U^ p -u n a ’ np \ 6 ™ less or equal than 2\a\h(E( Q n K \r%)+ 
E{g'i / \r r l)), and than 2\a\h(\U% K — u\ | 6 + \Ujf L — u^\ 6 ), respectively. Consequently,we get by the same 
reasoning as in the previous step, under assumption 7 > 3/2, 


\V 2 \<c(6 ,M 0 ,E 0 ,r,rA\(U,VU)\\ L ^ {QT . Rl2) )kY i £(Q n ,u n \r n ,U n ) + 5kY / \u n -Uf:)\U (m3y (8.4) 


n —1 


n =1 
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Step 3: Term V 3 . Since the pair ( r,U ) satisfies continuity equation (11 . laf) in the classical sense, we 
have for all n = 

[d t r] n + U n ■ Vr n = -r n divt/ n , 

where we recall that [3 t r] n (x) = dtr(t n ,x) in accordance with (13.171) . Using this identity we write 

m 

^ = ^3,1+^3,2, V Z ,i = kY,n,i> 

n= 1 

with ^ 3,1 = - 5Z /. (p(&) ~ P( r K)(e?< - Tk) - p(Tk)) div U n dx 
KeT Jl< 

and 7^ 2 = £ [ / r * ~ 0 « p'(r n K )(u n - U n ) • Vr" dx. 

KeT Jr r l< 

Now, we apply Lemma 17.21 in combination with assumption (11,61) to deduce 

m 

\V 3 ,i\ <c\\divU\\ LO o {QT) kY / £(Q n ,u n \r n ,U n ). (8.5) 

n= 1 


Finally, the same reasoning as in Step 2 leads to the estimate 
1^3,21 < hc{M 0 ,E 0 ,r,r, |p , |c(|r 1 r])ll(Vr ) VL0lb~(fi ; R 9 )) 

m m 

+ c(5 , || r, r, b / |c , ([r : ,f])l|Vr|| L cx. ( f2. R 3 ) ) k £ £(*>", u n |r", U n ) PSk^ 

n =1 n=l 


IX 


- u r 


u 2 


(Q;R 3 ) ‘ 


( 8 . 6 ) 


Gathering the formulae (18.11) and (18.31) - (18.61) with 5 sufficiently small (with respect to /T). we conclude 
the proof of Lemma 18.11 □ 

Finally, Lemma 18. II in combination with the bound (14.131) yields 

m— 1 

+ ck J2 £(Q n ,u n \r n ,U n y, 

n= 1 

whence Theorem O is a direct consequence of the standard discrete version of Gronwall’s lemma. 
Theorem EU is thus proved. 


w 


£/ m ) < c h A + Vk + k + £(g°, 


ii°|r°, 


£/ u 


9 - Appendix: Fundamental auxiliary lemmas and estimates 

In this section we report several results related to the properties of the Sobolev spaces on tetrahedra and 
of the Crouzeix-Raviart (C-R) space. We refer to the book Brezzi, Fortin [3] for the general introduction 
to the subject. 

We start with the inequalities that can be obtained by rescaling from the standard inequalities on a 
reference tetrahedron of size equivalent to one. 

Lemma 9.1 ( Poincare, Sobolev and interpolation inequalities on tetrahedra). Let 1 < p < 00 . Let 
do > 0 and T be a triangulation of Li such that d > do where 6 is defined in h‘3.1\) . Then we have: 

(1) Poincare type inequalities on tetrahedra 

Let 1 < p < 00 . There exists c = c(do,p) > 0 such that for all K € T and for all v € W 1,P (K ) we 
have 

lb - v k \\lp(k) < ch\\Vv\\ LP ( K ), (9.1) 

Vcr € S(K), ||n - v a \\ LP ( K ) < ch\\Vv\\ LP (K)- (9.2) 
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(2) Sobolev type inequalities on tetrahedra 

Let 1 < p < d. There exists c = c(6o,p) > 0 such that for all K E T and for all v E W 1,P (K) we 


have 


where p* = ■jf—. 


\\v ~ v K \\ L v*( K ) < c\\Vv\\ LP(K) , 

Vo - € £{K), ||u - v a \\ LP * {K) < c||Vu|| iP(A -), 


(9.3) 

(9.4) 


(3) Interpolation inequalities on the tetrahedra 

Let 1 < p < d and p < q < p*. There exists c = c(9o,p) > 0 such that for all K E T and 
v E W 1,P (K) we have 

\\v-v K \\ li(k) < ch : \\Vv\\ LP(A';R d )J (9.5) 

\\V - V<t\\li(K) < chP \N v \\ L P(K-R d )^ ( 9 - 6 ) 


where \ = £ + 

Combining estimates (19.1119.61) with the algebraic inequality 




1/9 


(9.7) 


i=l 


i =1 


for all (a = (ai,..., aif) E R L , 1 < q < p < oo, we obtain the following corollaries. 


Corollary 9.1 (Poincare and Sobolev type inequalities on the Sobolev spaces). Under the assumptions 
of Lemma \9.1\. we have: 

(1) Poincare type inequalities on the domain 12 

Let 1 < p < oo. There exists c = c(Qo,p) > 0 such that for all v E PP 1,p (f2) we have 


v - ^IIlp(O) = ( 

:ei 

b v k\\ I Lp(K)) 

1/p 

— C/l | | Vw | | J^p (fJ;Rd) ) 

(9.8) 


AeT 





E 

ll' y — V <* Ilp(A') 

\ 1/P 

) < ch Vw iP (n ;8 d) 

(9.9) 


KeT o-ef(A') 

where v and v a are defined by CSD and (13.81) . 

(2) Sobolev type inequalities on the domain 12 

Let 1 < p < d. There exists c = c(6 o,p) > 0 such that for all v E VP 1,P (12) we have 

ll' u — ^IIap*(Q) < c||Vu|| LP(a) , 

(E E W V ~ V °W P LP* (K)) /P ^ C H Vv llLP(n;M d )- 

KG'T a££(K) 


(9.10) 

(9.11) 


(3) Interpolation inequalities on the domain 12 

Let 1 < p < d and p < q < p*. There exists c = c(9o,p) > 0 such that for all v E VP 1,P (12) we have 

l|v - ®IU«(n) < c/»^||Vu||iP (n)> ( 9 . 12 ) 

(E E W v ~ v A q Lq{K )) q <chP\\Vv\\ LP(n; R d )’ ( 9 . 13 ) 

Kd'T ae£(K) 

where - = - + . 

q p p* 
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Corollary 9.2 (Poincare and Sobolev type inequalities on 14). Under assumptions of Lemma \ 9. 11 there 
holds: 

(1) Poincare type inequality in 14(0): Let 1 < p < oo. There exists c = c(Qo,p) such that for all v E Vh, 

lb - ®ll lp(Q) < ch\v\ v r (n) , (9.14) 


(E £ 


V - Vn 


I LP(K) 


< ch\v\yp m . 


(9.15) 


KeT <j€.£(K) 

(2) Sobolev type inequality in 14.(0): Let 1 < p < d. There exists c = c(9o,p) such that for all v E 14(0), 

lb - ^IIlp* (n) < c blv?(n)> ( 9 - 16 ) 


E E 


v - v„ 


I LP*(K) 


) P * < c\v\v*(ny 


(9.17) 


K&T a££(K) 

(3) Interpolation type inequalities in 14(0) 

Let 1 < p < d and p < q < p*. There exists c = c(0Q,p) > 0 such that for all v E 14(0) we have 


E E 

KgT a g£(K) 


\v — V. 


L1(K)J — 


W)’ 

(9.18) 

< ch^\v\ V P {n) , 

(9.19) 


where 1 = £. + LA' 


The next fundamental lemma deals with the properties of the projection Vh defined by (13.91) . 

Lemma 9.2 (Projection on 14)- Let 9$ > 0 and T be a triangulation of 0 such that 9 > 9 q where 9 is 
defined in (Eiy. 

(1) Approximation estimates on the tetrahedra 

Let 1 < p < oo. There exists c = c(9o,p) > 0 such that 

w e wt’ p (n)nL°°(n),\/K e t, \\v h \\ L °o {K) < c\\v\\ Lao{K) , 


Vu E wZ*nW a *(n),VK E T, \\v - v h \\ LP{K) < ch s \\V s v\\ LP{K . Rdy 
||V(u — v h)\\LP(K-,M. d ) — c h 11^ ^ 11 Lp (A’ ;R d< *) > ® = ^ ^' 

(2) Preservation of divergence 

Vv E W 0 1 ’ 2 (O,M d ),Vg E Lft(O), [ q diva^dx= [ q div-udx 

Jn 

(3) Approximation estimates of the Poincare type on the whole domain 

Let 1 < p < oo. There exists c = c(6o,p) > 0 such that for all v E Wq’ p (0), 


(9.20) 

(9.21) 

(9.22) 

(9.23) 


b-Oil|ip(n) < c/i||Vr:|| LP(n . ]f 


(9.24) 
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(4) Approximation estimates of the Sobolev type on the whole domain 

Let 1 < p < d. There exists c = c(6o,p) > 0 such that for all v E Wq’ p (Q), 

11^ — v h\\LP*(Q.) — c ll^ l 'llLP(n;R <i )- (9.25) 

Statement (2) of Lemma 19.21 is proved in |6j, where one can find also the proof of item (1) for p = 2. 
We present here the proof of statements (1), (3), (4) for arbitrary p for the reader’s convenience, since 
a straightforward reference is not available. 


Proof. Step 1: We start with some generalities. First we complete the Crouzeix-Raviart basis (13.61) by 
functions f> a indexed also with a E £ ex t saying 

S = S a y, (■ <J,a')e£ 2 

|cr | J a > 

and observe that 

^ 4>a{x) = 1 for any x E K. (9.26) 

a&£(K) 

A scaling argument yields 


|| 0 O-||x,°°(f2) — c (^o)> ^'ll'^ < ^0-|lL oo (0;R‘ i ) — c (^o)- 
Second, we define the projection v —>• Vh for any v E W 1,p (Ll) by saying 

Vh = y v a (j) a . 

cr££ 


(9.27) 


We notice that if v E Wq’ p (Q) then Vh coincides with (13.51) . Moreover, 

Vh = v for any affine function v. 


(9.28) 


Third, due to the density argument, it is enough to show the remaining statements (1), (3), (4) for 
v E Wq’ p (LI) fl W s,00 (n), s = 1,2, according to the case. 


Step 2: We realize that suppti^ = K U L and derive (19.201) directly by employing representation 
(ES|), definition of Vg- and estimate (19.271) . 

We denote by xk = f K xdx the center of gravity of the tetrahedron K. We calculate by using 
(|9.28l) and the first order Taylor formula 


v(x) - v h (x) = v(x) - v(x K ) - [v - v(x K )]h(x) 


= (x-Xk)- Vv(xk +t(x — XK))dt— y (f>a(x)- 
Q (Tg £(K) 


(J\ 


(x-xk)- / Vv(xk + t(x — XK))dt dS, 


where x E K. This formula yields immediately the upper bound stated in (19.211) ,_i if p = oo. If 
1 < p < oo we calculate the upper bound of the L p -norm of each term at the right-hand side separately 
by using (19.271) . Fubini’s theorem, Holder’s inequality and the change of variables y = xk + t(x — xk) 
together with the convexity of K. 

The same reasoning can be applied to prove (19.211) o-o. Indeed, we observe that 


v(x) - Vh(x) = v(x) - (x - x K ) ■ Vu(x^) - v(x K ) - [v - {x - x K ) ■ Vv(x K ) - v(x K )]h(x) 
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by virtue of (19.281) . Now we apply to the right hand side of the last expression the second order Taylor 
formula in the integral form, and proceed exactly as described before. 

Finally, one applies the same straightforward argumentation to get (19.221) . This completes the proof 
of statement ( 1 ). 

Step 3: Statement (3) follows easily from (I9.21lb -i and the algebraic inequality (19.71) . 

Step 4: We use (I9.2(il) and (19.281) to write 

v(x) - v h (x) = ^2 (v(x) - v a )</> a (x), x e K; 

aGS(K) 

whence 

11 ^ - v h\\LP* (K) — LP(K;M. d ) 

where we have used the Sobolev inequality (19.41) on the tetrahedron K E 7~ and the L°°-bound (19.271) . We 
conclude the proof of statement (4) by using the relation (19.71) . The proof of Lemma l9.2l is complete. □ 

The following corollary is a direct consequence of (19.221) . 

Corollary 9.3 (Continuity of the projection onto 14). Under assumptions of Lemma 1 9.21 there exists 
c = c( 6 o,p ) > 0 such that 

Vv E Wq’^H), b/i|yP(Q) 4 c ll^' u ll,LP(a;R d )’ (9.29) 

where 1 < p < oo. 

Although the non conforming finite element space 14 is not a subspace of any Sobolev space, its 
elements enjoy the Sobolev type inequalities. This important fact is formulated in the next lemma. 

Lemma 9.3 (Sobolev inequality on 14). Let Ut be a bounded domain ofM d . Let T be a triangulation of 
the domain f l in simplices such that 6 > 6 q > 0 where 0 is defined in (ED- Then we have: 

(1) Sobolev inequality in 14(H) (case 1 < p < d): 

There exists c = c(0q,p) such that for all v E 14(H), 

IMIlp*( 0 ) ^ c Mvjf(n)- ( 9 -30) 


(2) Sobolev inequality in 14(H), case p > d 

Let 1 < q < oo. There here exits c = c(9o,p,q) > 0 such that forall v € 14(H), 

IMU-qn) < c\v |v£(fi) (9-31) 


Proof. Step 1 Let 1 < r < a < oo. Let u E 14- We call v the element of 14 such that v a = \u a \ a . Then 
there exists C only depending on d, r, a such that 


\u\ 


(«) - 


\u r r 


L*(ny 


(9.32) 


To prove (19.321) we remark that, using a change of variable, it is enough to show to prove the existence 
of C for only the unit symplex K. Let u E Pi (A') and we call v the element of Pi (A') such that v a = |ix CT | a . 


Let Tiu) = I \u\ \r r (K\ an d S(u) = ||«|| “r ~ . These two functions are continuous, homogeneous of degree 
^ V 2 *-) L ol (K) 

1 and non zero if u 0. Since Pi (AT) is a finite dimensional space, we can choose a norm on Pi (A") and 
take C = ( 7 ) 7 )“ where M = max{T(ii), ||rt|| Pi ^ = 1} and m = min{T(it), ||rt|| p /^ = 1}. 
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Step 2: Proof for p = 1. 

We set u = 0 outside fL For cr E £j nt , cr = Ji|L, we set |[u(x)]| = |u^(x) — ul(x)\ for x E cr. For 
cr E 5ext n£(iF), we set |[rt(x)]| = \uk(x)\ for x € cr. We first remark that there exists Cgi and C 1.2 only 
depending on d such that 


Ml d < Cl ll 

Ld-l(O) 


Mlw(R d ) < CiMIV/MIl 1 ^) + Ci,2 E f |M|d5- 

a£S Ja 


We now prove that there exits Cg 3 only depending on d and 6 q such that 

E / IM|d5<C' 1 , 3 ||V fc t i || L i (n) . 


Let K E T and cr E £{K). Let x a be the center of mass of cr. We have, wittm/< = u in K, 

ur{x) — u(x a ) = / Vuk ■ (x — x a ) dx. 

Jo 

Then if cr = K\L we have 

\uk(x) - u L (x )| < h a (\Vu K \ + |Vul|). 

Integrating this inequality on cr gives 

J IMI d5, < \cr\h c (\Vu K \ + |Vu L |) < ^(||Vu|| L i w + ||Vu|M ( l) ). 

Similarly for cr E £ ex t Fl £{K) we have 

J IMId>5 < |Vn| \ L i( K ) 

Then there exists Cl^ = C(d, # 0 ) such that 

E / IM|dS , <C li 3||V ft «|M ( n). 

i,ae£ Ja 


and then, 

ll“llL 1 *(n) ^ c (^ 6 'o)||V /l n|| L i (0) . 

Step 3: Proof for 1 < p < d. 

Let 1 < p < d and p* = and let u E VC We set u = 0 outside fb Let a = p ^_.^ , so that a > 1 and 
al* = p*. We call v the element of 14 such that v a = |wo-| a for cr E £. One has v / |n| a but there exits 
C 24 only depending on d and p (see lemma IU. 321) such that 


IMI2p* ( n) ^ c 2 I i||v||z,i* ( n) < c(d,p,0o)||V ft u|| L i (n) . 
Moreover using a scalling argument we obtain 

\\Vhv\\mK)<c(d, P ,e 0 ) E l^r-'lv^lMl. 

ae£(K) 

Then, using Holder Inequality, we have, with q = (so that q(a — 1) = p*), 

P* 

l|Vftu||z,i ( if) < c(d, p, 0 O )| | Vn| \lp(k) 1M l/ P * (ft:) • 
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Summing on K £ T we obtain 


w IIlp*(H) ^ C 2 WV h u\\ L p^ n y 


Step 4: Proof for p > d. 

Let 1 < q < 00 . There exists r = r(d , q) such that r < d and r* > q. We have 

IMIz/-*(n) < c(r,d,q,9 0 )\\V h u\\ L r {n y 

Moreover 


MU«(n)<|H|« ^IMIz/Rn) ^ c(d,q,9 0 )\n\i '■*||V h u|| L r (n) 


and 

Finally 


IIVHIl^q) < |^| r 

u\\l<ksi) < c(ft,d,p,q,eo)\\V h u\\ LP ( n ). 


□ 


A Combination of Lemma 19.31 with estimates (|9.14|) . (19.161) and the Holder inequality yields the 
following corollary. 

Corollary 9.4 (Estimates of the norms of mean values). We have under the assumptions of Lemma \ 9. dl 

(1) Poincare type inequality involving mean values on tetrahedra 

There exists c = c(9q,p ) such that for all v € 14, 

ll®|| LP(SI) = ( \ K W v k\ p ) /P < c(|M| L r(fi) + h\v\ V P {n) ). (9.33) 

K&T 

(2) Sobolev type inequality involving mean values on tetrahedrons 
Let 1 < p < d, there exists c = c(9$,p) such that for all v € 14, 

II^IIlp*( p) = ( II \K\\vk\ p *) /P < c(||u|| LP . (n) + Mvjp). (9.34) 

K&T 


Note that the Last but not least, we recall a result on equivalence of norms in the space 14(f)) which 
is a consequence of a discrete Poincare inequality on the broken Sobolev space 14 [331 proposition 4.13]. 

Lemma 9.4 (Discrete and continuous norms in 14). Let 1 < p < 00 . Let 9 q > 0 and T be a triangulation 
of Cl such that 9 >9 q where 9 is defined in (ED- Then the norms 

( \ a \ h \ v v\ P ) /P and IMIlp(0 ) (9.35) 

0"(=£int 

are equivalent on 14 (H) uniformly with respect to h > 0. 


The last lemma in this overview deals with the estimates of jumps over faces. The reader can 
consult [8, Lemma 3.32] or [IS, Lemma 2.2] for its proof. 


Lemma 9.5 (Jumps over faces in the Crouzeix-Raviart space). Let 9o > 0 and T be a triangulation of Cl 
such that 9 > 9 q where 9 is defined in ED- Then there exists c = c{6q) > 0 such that for all v € 14(H), 



(9.36) 


where [v\ a , na is a jump of v with respect to a normal n a to the face a, 


and 


\Zx £ (7 - K\L £ £intl It] (7, Tier ( 3 O 


v\k(x) - v\l(x) if n a = rifj K 
v\l(x) - v\k(x) if n a = n^ L 


Vx € a £ £ ex t) [ v \<r,n a (x) = v(x), with n a an exterior normal to dCl. 
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